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Introduction 


For  improved  care  of  patients  with  acute  lung  injury,  existing  clinical  data  sets  are 
examined  to  characterize  quantitatively  the  pulmonary  pressure- volume  (p-V)  curves.  A 
mechanistic  model  of  the  total  respiratory  system  (TRS)  is  constructed  for  the  inflation 
process,  making  it  possible  to  relate  the  p-V  curve  to  the  corresponding  changes  in  intra- 
respiratory  conditions.  The  mechanistic  model  computationally  simulates  the  TRS  based 
on  an  application  of  the  principle  of  statistical  mechanics  to  a  very  large  nmnber  of  ele¬ 
ments  comprising  the  TRS.  A  mechanistic  model  of  the  deflation  process  is  currently  imder 
investigation. 


Body 


Research  contents  are  reported  in  three  chapters: 

Chapter  1.  Quasi-Static  Pressure- Volume  Curce:  Comprehensive  Data  Analysis, 

reports  research  results  corresponding  to  Objective  1  (Examination  of  accuracies 
and  limitations  of  the  sigmoidal  equation)  and  Objective  2  (Development  of  a 
method  for  quantitative  characterization  of  p-V  curves)  in  Statement  of  Work  of 
Research  Proposal. 

Chapter  2.  Mechanistic  Model:  Part  I.  Model  Development  for  Inflation  Process, 
reports  the  derivation  of  the  mechanistic  model  of  the  total  respiratory  system, 
which  is  part  of  Objective  3  (Development  of  a  mechanistic  respiratory  model)  in 
Statement  of  Work  of  Research  Proposal. 

Chapter  3.  Mechanistic  Model:  Part  II.  Examination  of  Clinical  Data, 

reports  examinations  of  p-V  inflation  curves  based  on  the  mechanistic  model. 
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Chapter  1.  Quasi-Static  Pressure- Volume  Curve: 

Comprehensive  Data  Analysis 

Abstract 

A  p-V  model  equation  with  four  parameters  is  used  to,  represent  various  existing  (p-V) 
curves.  The  report  is  focused  on  the  case  in  which  the  equation  is  applied  to  two  existing 
groups  of  p-V  data  (one,  twenty  nine  p-V  cmrves  of  healthy  adults  and  the  other,  twenty 
one  p-V  curves  of  patients  with  acute  respiratory  distress  syndrome)  to  determine  the 
magnitudes  of  the  parameters  for  each  data  set.  The  equation  is  found  to  represent  the  p- 
V  curves  of  both  data  groups  extremely  well.  It  is  also  confirmed  that  the  magnitudes  of  the 
four  parameters  of  the  error  function  p-V  model  equation,  combined  with  the  corresponding 
normalized  representation  of  p-V  curves,  quantitatively  distinguish  different  respiratory 
conditions  between  the  two  groups  as  well  as  between  different  data  sets  in  each  group. 
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Nomenclature 


p 

pressure  (interpleural  pressure  difference) 

Pgrad 

(volume-)  gradient  pressure  range,  Eq.(3) 

Pcl{u) 

lower  (upper)  corner  pressure,  Eq.(4) 

Pmci{d) 

pressure  at  maximum  compliance  increase  (decrease),  Eq.(5) 

Po 

pressme  at  the  inflection  point  (at  the  maximum  local  compliance) 

where  V  =  {Vv  V£,)/2 

P 

non-dimensional  pressure,  p/po  —  1 

V 

volume 

Vl{u) 

lower  (upper)  volmne  asymptote  (Fig.l) 

V 

non-dimensional  volume,  (V  —  ({Vu  -I-  Vjr,)/2)/(AV/2),  (Eq.(lb)) 

AF 

Vu  -  Vl  (Fig.l) 

Greek  symbols: 

a 

constant  deflned  in  Eq.(la) 

A 

apo A V’(non-dimensional  parameter)  (Eq.(lb)) 

0) 

Ap/2  (Eq.(lb)) 

Acronyms: 

ARDS 

acute  respiratory  distress  syndrome 

LIP 

lower  inflection  point,  Eq.(6) 

TRS 

total  respiratory  system 

UIP 

upper  inflection  point,  Eq.(6) 
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Introduction 

Quasi-static  pulmonary  p—V  (pressure  -  volume)  curves  provide  quantitative  informa¬ 
tion  on  the  respiratory  system  that  is  important  for  both  research  and  clinical  guidances. 
A  typical  inflation  p-V  curve,  obtained  for  an  anesthetized  human  subject  in  supine  posi¬ 
tion,  consists  of  a  nearly  linear  region  of  high  compliance  (i.e.  large  dV/dp)  sandwiched 
between  two  segments  with  low  comphance  at  low  and  high  pressure  regions.  The  shape 
of  the  curve  is  affected  by  two  mechanisms,  the  distension  of  the  elastic  respiratory  wall 
tissue  components  and  the  recruitment  of  the  alveoh  (‘pop-open’  mechanism).  The  latter 
is  the  opening  of  alveoU  overcoming  the  surface  tension  at  the  interface  between  the  gas 
and  the  liquid  film  lining  the  alveolar  surface.  A  pressure  increase  (i.e.  an  increase  in  the 
interpleural  pressure  difference)  results  in  the  recruitment  of  a  greater  number  of  alveoli. 
The  high  compliance  is  believed  to  be  associated  with  both  the  distension  of  open  parts 
and  the  (alveolar)  recruitment  of  collapsed  parts  of  the  total  respiratory  system  (TRS). 

In  order  to  quantify  the  characteristics  of  p  —  V  curves  as  well  as  their  changes  ob¬ 
served  in  clinical  settings,  various  p-V  model  equations  have  been  proposed  [1-8].  One 
commonly  used  model  equation  is  developed  by  dividing  the  entire  p-V  curve  into  three 
regions,  a  high-pressmre,  low-compliance  upper  region,  a  high-compliance  midregion  and 
a  low-pressure,  low-compliance  lower  region.  The  midregion  is  represented  by  a  linear 
equation  between  p  and  V;  while,  the  two  low-compliance  regions  are  approximated  by 
an  exponential  function  of  pressure  [9,  10].  The  linear-exponential  model  equation  is  a 
piecewise  continuous  function  with  the  compliance  abruptly  changing  its  magnitude  at  the 
intersects  of  the  linear  and  the  exponential  regions.  Venegas,  Harris  and  Simon  [8],  on 
the  other  hand,  showed  that  a  single  continuous  function  in  a  form  of  sigmoidal  (tangent 
hyperbolic)  equation  represents  various  p-V  curves  extremely  well.  Parameters  in  model 
equations  (both  piecewise-continuous  and  continuous  equations)  are  determined  from  sta¬ 
tistical  processing  of  clinical  data.  More  recently  the  clinical  usefulness  of  the  sigmoidal 
model  equation  over  piecewise-continuous  representations  is  also  reported  by  the  same 
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group  of  researchers  [11]. 

Accurate  and  quantitative  determinations  of  the  form  of  p-V  model  equation  and 
its  parameters  are  prerequisite  to  the  clinical  interpretations  of  p-V  curves,  including  an 
establishment  of  ventilator  strategy  with  the  p-V  curve  guidance  in  intensive  care  for 
patients  with  acute  lung  injury  as  well  as  its  more  severe  form,  acute  respiratory  distress 
syndrome  (ARDS)  [12  -  14].  Om  objective  is  to  test  a  hypothesis  that  the  continuous 
p-V  model  equation,  particularly  in  a  form  of  an  error  function  equation,  is  effective  in 
representing  p-V  cmves  from  different  sovuces,  and  of  different  respiratory  conditions  (a 
group  of  patients  with  ARDS  [8]  and  a  group  of  healthy  adults  [9]).  The  former  covers  both 
inflation  and  deflation  processes,  and  the  latter  includes  the  inflation  p-V  curves  before  and 
after  alveolar  recruitment  maneuver  with  a  total  of  fifty  p-V  curves.  The  report  examines 
differences  and  similarities  (1)  between  patients  with  ARDS  and  healthy  adults,  (2)  among 
patients  with  ARDS  as  well  as  among  healthy  adults,  both  in  terms  of  parameters  of  the 
error  function  p-V  model  equation. 

Equations  for  Quasi-Static  p-V  Curves 

A  model  equation,  originally  proposed  by  Venegas,  Harris'  and  Simon  [8]  and  subse¬ 
quently  shown  to  represent  p-V  curves  well  for  both  illation  and  deflation  processes  [11], 
has  the  following  sigmoidal  (tangent  hyperbolic)  form; 

^  =  -a(V-Vu)(V-VL),  ^^^^==[l  +  exp{-a/^V{p-po))]-\  (la) 


where  AV  =  Vu  —  Vi,,Vu  —  upper  volume  asymptote,  Vl  =  lower  volume  asymptote,  a  = 
positive  constant  and  po  =  pressiue  at  the  midpoint  (inflection  point)  of  the  curve.  The 
corresponding  non-dimensional  form  of  the  sigmoidal  equation  is  [15], 

—  =  -  — (V.-l),  V  — - tanh{uj))  (16) 


where 


V  -  ^  -  (^t;  +  Vl)/2 
AV/2 


P  ^ 

Po 


A  =  apo  AV. 
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Venegas,  Harris  and  Simon  suggested  that  a  p-V  equation  in  terms  of  the  error  function 
is  also  effective  in  representing  p-V  data  [8]. 

The  error  function,  erf(x),  is  defined  as 

er/(a:)  =  — 7=  /  e~*  dt  with  erf{oo)  =  l,  erf{—x)  =  -erf{x). 

v’T  Jo 

The  error  function  model  equation  may  be  expressed  as 

%  =  ^^exp[-(^aAK(p  -  p„))'],  =  i[l  +  er/(^Ap)] ,  (2a) 

^  =  exp{-  ^  u;^) ,  V"  =  erf(^w) .  (26) 

Fig.l  is  a  sketch  of  a  typical  p-V  model  equation  (either  the  sigmoidal  or  the  error-function 
model  equation).  The  cmrve  varies  smoothly  between  the  low  pressure  asymptote,  V/,,  and 
the  high  pressure  asymptote,  Vu.  The  midpoint  of  the  curve  where  the  volume  is  equal  to 
(Vu  +Vl)/2  is  the  inflection  point  with  its  pressure  denoted  by  po.  Both  the  sigmoidal  and 
the  error-function  model  equations  are  antisymmetric  with  respect  to  the  inflection  point; 
that  is,  V{p  -  Pq) -V{po)  = -{V(po  -  p)  -  V{po))  or  F(p)  =  -V{-p).  The  compliance, 
dV/dp,  increases  along  the  p-V  equation  as  pressure  increases,  until  the  inflection  point 
(=  the  point  of  maximum  compliance)  is  reached.  Then  the  compliance  decreases  with  a 
further  increase  in  pressure.  A  tangent  to  the  model  equation  curve  at  the  inflection  point 
has  the  compUance  of  Q!(AV)^/4.  The  two  points  of  intersection  between  the  tangent 
and  the  two  volume  asymptotes,  V  —  Yu  and  V  =  Vl  axe  referred  to  as  the  upper  and 
lower  corner  pressure,  Pcu(cZ)>  respectively.  The  pressiure  difference  between  the  two  corner 
pressures  is  defined  as  the  (volume-)  gradient  pressure  range,  Pgrad-  Also,  the  pressure  at 
the  point  of  maximum  compliance  increase  (decrease)  of  the  p-V  curve,  Pmd  (Pmcd ),  may 
be  specified  as  the  points  where  the  third  derivative  of  V  with  respect  to  p  is  zero. 

For  both  the  sigmoidal  and  the  error-function  model  equations. 


Pgrad  f  _  _ ^  ^ 

Po  y  Po{dV/dp)  max.  J 


Pcu(cl) 

Po 
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(3,4) 


On  the  other  hand, 

_  ! _  Pmcd^mci)  ,  _  J  (~)1- 

Pmcd{mci)K—  ~  ~ 

Fig.2  is  the  (non-dimensional)  p  —  V  curve,  corresponding  to  the  p-V  curve  of  Fig.l. 
The  origin  (p  =  0,  F  =  0)  represents  the  point  of  antisymmetry  {po,  {Vu  +  Vi)/2)  of  Fig.l. 
The  non-dimensional  pressure,  p,  is  the  pressure  difference,  p  -  po,  as  a.  fraction  of  po 
(Eq.(lb).  The  normalization  of  volume  shifts  the  upper  and  the  lower  volume  asymptotes, 
Vu  and  Vl  in  Fig.l  into  -f-1  and  —  1  respectively  in  Fig.2.  With  both  the  location  of  po 
and  the  volume  asymptotes  made  common  to  all  p-V  curves,  the  resulting  non-dimensional 
representations  characterize  p)-V  curves  in  general  in  terms  of  a  single  non-dimensional 
parameter,  A.  (Eqs.(lb,2b))  From  Eq.(3)  the  parameter.  A,  is  four  times  the  ratio  of  the 
pressmre  at  the  maximtun  compliance,  po,  to  the  volume-gradient  pressure  range,  Pgrad- 
Since  the  compliance  is  maximum  at  the  origin,  the  first  quadrant  (V,  p  >  0)  in  Fig.2  is  a 
region  of  decreasing  local  compliance  with  pressure;  while,  the  third  quadrant  (V,  p  <  0) 
is  a  region  of  increasing  local  compliance  with  pressure.  The  origin  (p  =  0,  V  =  0)  of 
dimensional  p-V  curves  is  transformed  into  {p  =  —1,  V(V  =  0))  on  a  p-V  curve;  hence, 
the  physiological  lower  limit  of  p  is  —1.  Various  pressure  locations  on  p  —  V  diagram  are 
proportional  to  1/A  as  shown  in  Eqs.(3-5).  Eqs.(3-5)  also  imply  over  the  pressxire  range 
of  p  >  0  that  there  is  no  lower  corner  pressure  (i.e  1  +Pci  <  0)  if  A  <  2,  and  that  there 
is  no  pressure  for  maximum  compliance  increase  (i.e.  1  -f-p^ci  <  0)  if  A  <  1.317  (1.596) 
for  the  sigmoidal  (error-function)  model  equation.  Both  the  sigmoidal  and  the  error- 
function  model  equations  are  capable  of  representing  p-V  curves  over  their  entire  ranges 
as  continuous  functions. 

Piecewise-continuous  model  equations  are  also  used  to  represent  p-V  curves.  Shown 
below  is  a  three-region  model  equation  [9],  relevant  to  the  present  study,  consisting  of  a 
linear  midregion  (Vl/p  <  V  <  Vuip)  and  two  exponential  regions  at  high  (Vy/p  <  V)  and 
low  pressure  (V  <  Vup)  ranges.  (The  subscripts,  L(U)IP  (lower  (upper)  inflection  point). 


3 1 7/ A  for  sigmoidal  equation; 
596/ A  for  error-function  equation. 
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indicate  the  points  where  the  linear  midregion  equation  intersect  with  the  exponential 
equations.) 

Vmax  -  V (p)  -  (Vmax  “  V[//p)  •  (p  -  puip)^ 

for  VuiP  <v  <  Vinsuff,  - 
V{p)  -  Vlip  =  Clin  (p-PLip)  for  Vlip  <V  <  Vuip,  (6) 

V„i„-V(p)  =  (f-pup))  for  V  <Vup, 

where  Cun  =  compliance  at  the  hnear  midregion,  Vmin^max)  =  volume  asymptote  of  the 
lower  (upper)  exponential  region, 

Data  Analyses 

The  two  data  sources  with  a  total  of  50  p-V  data  sets  consists  of  (A)  ARDS  patients 
by  Harris  et  al  (2000)  [11],  and  (B)  healthy  adults  by  Svantesson  et  al  (1998)  [9].  For  the 
data  source  A,  p-V  data  points  were  made  available  to  us  by  the  authors.  The  data  source 
B  provides  model  parameters  of  the  piecewise-continuous  model  equation,  Eq.(6),  as  well 
as  data  ranges  for  each  data  set.  Information  on  the  data  soiurces  relevant  to  the  present 
study  is  summarized  below. 

Data  Source  A 

21  data  sets  of  ARDS  patients  by  Harris,  Hess,  Venegas  [11],  Original  p-V  data  points 
made  available  by  the  authors.  Inflation  and  deflation  data  in  supine  position. 

Data  Somce  B 

29  data  sets  of  healthy  adults  (both  male  and  female)  by  Svantesson,  Sigurdsson, 
Laxsson,  Jonson  [9],  14  data  sets  before  and  15  data  sets  after  alveolar-  recruitment 
maneuver.  Inflation  data  in  supine  position 

The  parameters  of  Eq.(6)  (Vmin,  Vmax,  Vlip,  Vuip,  Vnauff,  Cun,  Pup,  Puip)  are 
tabulated  for  all  data  sets  in  [9] . 

Data  sets  from  the  somrce  A  axe  analyzed  by  minimizing  the  difference  between  data 
points  and  the  model  equation  (either  the  sigmoidal  or  the^error-function  model  equation) 


12 


through  the  application  of  the  method  of  least  squares  to  obtain  the  parameters,  A,  A  V,  po 
and  Vu  (or  Vl).  To  analyze  the  data  source  B,  ten  to  twenty  five  computational  data 
points,  depending  on  the  data  range,  are  generated  from  Eq.(6).  Then,  the  method  of 
least  squares  is  applied  to  determine  the  parameters.  A,  poi  AV  and  Vu  (or  Vl)  of  the 
error  function  model  equation.  Parameters  of  p-V  model  equations  are  determined  for  all 
data  sets  in  Data  Source  A  (for  both  inflation  and  deflation  data)  and  B  (for  both  before- 
and  after-recruitment  maneuver).  Discussion  beyond  the  validity  test  of  the  error  function 
model  equation,  however,  is  focused  on  the  inflation  data  sets  of  the  two  data  somrces. 
Results  of  data  analyses  are  summarized  in  Table.  1  (for  Data  Source  A)  and  in  Table.2 
(for  Data  Source  B). 

Results  and  Discussion 

Fig.3  shows  a  typical  data  set  of  an  ARDS  patient  from  the  data  source  A  as  well 
as  the  sigmoidal  and  the  error  function  model  equations,  Eqs.(la)(2a),  determined  by  the 
method  of  least  squares.  The  parameters  of  the  model  equations,  (A,  po)  Vu),  are 
(1.470,  13.308,  3.491,  2.750)  for  the  sigmoidal  equation  and  (1.627,  13.324,  3.156,  2.584) 
for  the  error  function  equation.  Both  equations  represent  the  data  points  well  over  the 
entire  data  range.  Substantial  differences  between  the  two  equations  occur  in  high  and 
low  pressure  regions  away  from  the  data  range  as  they  approach  different  asymptotes  of 
V  =  Vu  and  V^.  It  should  also  be  noted  that  there  is  no  lower  corner  pressure  for  the  data 
set  since  A  is  less  than  2  (Eq.(4)),  and  that  the  pressure  at  maximum  compliance  increase, 
Pmcii  is  very  low  at  1.385  cmH20  for  the  sigmoidal  model  equation  and  at  0.254  cmH20 
for  the  error  function  equation  (Eq.(5)). 

An  example  of  the  analysis  of  the  data  source  B  is  shown  in  Fig.4.  The  dotted  curve 
represents  the  piecewise  continuous  equation,  Eq.  (6),  with 

{Vmin,  Vmax,  VliP,  VuiP,  Vinsuff,  Clin,  PLIP.i  PUIp) 

=  (-2230  [mL],  5870,  1513,  2884,  4125,  157  [mL/cmRzO],  14.4  [cmH^O],  23.1), 
reported  in  [9].  The  solid  curve  is  the  corresponding  error  function  model  equation,  Eq.(2a), 
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obtained  by  applying  the  method  of  least  squares  to  computational  data  points  generated 
over  the  data  range  of  0  <  V  <  Vinauff-  Since  “true”  data  points  are  unknown  in  the 
Unear  region  of  Eq.(6),  we  imposed  a  constraint  that  the  inflection  point,  po>  of  the  error 
function  equation  is  located  in  the  midregion  of  Eq.(6).  The  parameters  of  the  error 
function  equation  thus  determined  are  (A,  poi  ^u)  =(2.0800,18.224,5.58750,  4.9114). 
Due  to  the  error  minimization  Eq.(2a)  is  nearly  identical  to  Eq.(6)  over  the  data  range 
indicated  by  the  two  triangle  malrks;  which  is  valid  for  other  28  data  sets  from  Data  Source 
B.  A  continuous  change  of  the  compUance  (i.e.  non-linear  p-V  change)  in  the  region  near 
Po  has  been  described  previously  in  terms  of  the  sigmoidal  (tangent  hyperbolic)  model 
equation  [15].  It  should  also  be  mentioned  here  that,  of  twenty  nine  inflation  data  sets, 
Po  (the  inflection  point)  is  between  plip  and  puip  in  eighteen  data  sets,  equal  to  puip  in 
eight  data  sets,  equal  top  up  in  one  data  set,  and  po  =  PuiP  =  Plip  in  two  data  sets. 

Fig.5.  is  a  plot  of  the  error  function  equation,  Eq.(2a),  with  -^Tr  A(p/po  —  l)/4  and 
(V^  —  Vl)/AV  as  X-  and  y-  axis  respectively.  All  (both  inflation  and  deflation)  data  points 
of  Data  Source  A  are  also  shown  in  the  figinre,  confirming  very  good  agreements  with 
the  equation.  The  coefficient  of  determination,  B?,  is  0.999247,  which  is  comparable  in 
magnitude  to  that  of  the  sigmoidal  equation,  =  0.9992  reported  in  [11],  thus  indicating 
that  both  the  sigmoidal  (tangent  hyperbolic)  p-V  equation,  Eq.(la),  and  the  error  function 
equation,  Eq.(2a),  are  very  effective  in  representing  quasi-static  p-V  curves.  Shown  in  Fig.  6 
are  comparisons  between  the  sigmoidal  equation  and  the  error  function  equation  in  terms 
of  two  parameters  in  the  equations,  A  (Fig.6(a))  and  po  (Fig.6(b))  for  twenty  one  data  sets 
from  Data  Source  A.  Due  to  differences  in  functional  form  the  magnitude  of  A  is  slightly 
higher  for  the  error  function  equation  than  for  the  sigmoidal  equation.  On  the  other 
hand,  the  inflection  point,  po,  being  the  point  of  antisymmetry,  should  be  identical  in 
theory  for  both  model  equations.  Fig.6(b)  confirms  it  as  the  magnitudes  of  po  determined 
by  the  method  of  least  squares  are  very  close  between  the  two  equations.  Differences  in 
the  magnitudes  of  A  and  po  between  the  two  continuous-function  model  equations  would 
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result  in  function-specific  values  for  such  quantities  as  Pcu{ci)y  Vmcd(mci)  of  Eqs.(4,5)  which 
characterize  p-V  curves,  indicating  the  importance  of  using  the  same  p-V  model  equation 
in  order  to  analyze  clinical  data  in  a  consistent  maimer. 

Although  quantitative  comparisons  of  parameters  cannot  be  made  among  different 
p-V  model  equations  in  a  mathematically  rigorous  manner,  a  comparison  between  the 
Unear-exponential  model  equation,  Eq.(6),  and  the  error  function  equation,  Eq.(2a),  is 
presented  in  Fig.7  in  a  form  of  Pmd  vs  piip  in  Fig.7  (a)  and  Pmcd  vs  puip  in  Fig.7  (b)  for 
Data  Source  B.  It  may  be  seen  that  the  parameters  from  both  equations  distinguish  two 
data  groups,  Before  recruitment  maneuver  and  After  recruitment  maneuver,  successfully 
and  also  that  Pmcd  of  Eq.(2a)  and  puip  of  Eq.(6)  distinguish  the  two  groups  more  clearly 
than  Pmd  and  plip-  Values  of  Pup  {u IP)  must  be  located  directly  from  p-V  curves;  while, 
Pmd(mcd)  are  automatically  generated  from  the  model  equation  once  parameters  of  the 
model  equation  are  determined. 

As  shown  in  Eqs.(2a,b),  when  p  and  V  are  made  non-dimensional  the  resulting  non- 
dimensional  p-V  equation,  Eq.(2b),  contains  A  (=the  ratio  of  {pcu  —Pci)  to  po)  as  the  only 
parameter  representing  a  shape  of  p  —  V  curves  with  A/2  being  the  gradient,  dV /dp,  at 
the  origin  (where  p  =  0  i.e.  p  =  po.).  The  pdV  work  associated  with  the  process  from 
the  initial  to  the  end-of-infiation  pressure  was  suggested  as  a  quantity  representing  the 
pressure  range  actually  covered  by  a  specified  p-V  curve  [15].  However,  in  analyzing  data 
sets  from  different  sources  we  found  that  the  end-of-inflation  pressure  (volume)  is  data- (or 
investigator-)  dependent,  and  may  not  be  appropriate  as  a  comprehensive  indicator  for 
data  interpretation.  Here  we  selected  two  volume  differences  to  distinguish  p-V  data  sets 
accounting  for  the  range  of  data  relative  to  the  entire  range  covered  by  the  p-V  equation; 
one  is  Vb  —  V(p  =  0)  as  a  volume  scale  indicating  the  total  available  volume  of  a  specified 
TRS,  and  the  other  is  ^(p  =  20  cmH20)  -  V(p  =  0)  as  a  volume  scale  representing  volume 
range  covered  by  the  specified  TRS.  The  volume  at  p  =  20  cmH20  is  selected  arbitrarily; 
however,  in  all  the  data  sets  we  analyzed,  p-V  curves  were  obtained  beyond  p  =  20  cmH20. 
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Fig.8  shows  A  vs  po  in  (a)  and  [Vi;  —  V{p  =  0)]  vs  [V’(p  =  20  CTnH20)  —  V{p  —  0)]  in  (b) 
obtained  from  the  error  function  model  equation  for  all  inflation  data  sets.  In  terms  of  A, 
a  wide  range  (1  <  A  <  6)  is  covered  by  ARDS  patients;  while,  the  range  of  A  for  healthy 
adults  is  ~  1.5  —  3.5.  The  ^veolar  recruitment  maneuver  lowers  the  magnitude  of  A  of 
the  group  of  healthy  adults  as  a  whole.  A  similar  observation  may  be  made  on  po-  Both 
\yu  —  ^{P  =  0)]  and  [V{p  =  20  cmH^O)  —  V{p  —  0)]  in  Fig. 8(b)  are  low  in  magnitude  for 
patients  with  ARDS.  On  the  other  hand,  the  recruitment  maneuver  shifts  the  location  of 
the  whole  group  to  the  right  in  Fig.8(b).  Two  data  sets,  No.20  and  5,  representing  extreme 
points  in  Fig.8  (a)  and  (b)  respectively,  clearly  show  they  are  quite  different  from  those  of 
healthy  adults,  if  Fig.8(a)  and  (b)  are  examined  together. 

According  to  Data  Source  B  [9],  after  a  p-V  curve  before  the  recruitment  maneuver 
is  recorded,  the  lungs  are  inflated  to  an  airway  pressure  of  A0cmH2O  and  maintained  for 
15  s,  followed  by  six  pressure-controlled  breaths  (six  breaths/min.)  delivered  at  an  airway 
pressure  of  30  cmH20.  Then  a  second  large  insufflation  is  delivered  before  recording  a 
p-V  curve  after  recruitment  maneuver.  To  examine  the  p-V  curves  of  healthy  adults  as 
well  as  effects  of  the  alveolar  recruitment  maneuver  in  more  detail,  a  ratio  of  the  pressure 
at  the  inflection  point,  (po  (before  maneuver) /po  (after  maneuver)),  is  plotted  against 
a  ratio  of  A  {K{beiore  maneuver)  f  K  {Biter  maneuver))  in  Fig.9.  Each  data  point  is 
accompanied  by  two  numbers  indicating  the  data  set  number  and  his  or  her  age  in  the 
bracket  (unfilled  circle  for  male  and  filled  circle  for  female).  The  data  sets  of  the  younger 
may  be  seen  to  be  located  to  the  left  half  of  the  figure,  compared  to  the  older,  implying 
that,  for  the  healthy  young  adults,  A  after  the  maneuver  either  increases  slightly  or  remains 
roughly  the  same  as  A  before  the  maneuver.  In  order  to  discuss  implications  of  Fig.9 
further.  Figs.  10  and  11  show  the  error  ftmction  p-V  equations  before-  and  after-maneuver 
along  with  the  corresponding  non-dimensional  p—V  equations  for  three  data  sets  in  Fig.  10 
from  the  region  of  A  {before  maneuver)  /  h.  {after  maneuver)  <  ~  1.2,  and  for  three  data 
sets  in  Fig.ll  from  the  region  of  A  {before  maneuver) / h.  {after  maneuver)  >  ~  1.2.  The 
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equations  are  plotted  over  the  measurement  range  covered  by  Data  Source  B.  Numerical 
values  of  these  data  sets  are  tabulated  below. 


Data  No.  (Age)  1  (33)  6  (25)  7  (60)  11  (55)  13  (50)  15  (58) 

A  (before) /A  (after)  0.730  0.982  1.143  1.390  1.303  1.625 

Pq  (he f are) /pq  (after)  1.726  1.196  1.840  1.312  2.083  1.753 


Referring  to  Fig.  10,  the  high  po  ratios  of  Data  1  and  7,  compared  to  Data  6,  are  results  of 
substantial  reduction  in  po  after  the  maneuver  for  these  data  sets  as  observed  in  the  p-V 
equations.  The  triangidar  marks  on  the  p  —  V  curves  in  Figs.  10  and  11  indicate  locations 
of  p  =  2QcmH20',  hence,  onp  —  V  diagrams,  a  large  change  in  po  is  reflected  by  a  large 
shift  of  the  triangle  from  the  before-recruitment  location  to  the  after-recruitment  location. 
Different  degrees  of  changes  in  the  magnitude  of  A  over  the  alveolar  recruitment  maneuver 
for  the  three  data  sets  cannot  directly  be  observed  from  the  p-V  diagrams.  However,  on  the 
p-V  diagrams,  A/2  is  the  slope  of  p  —  V  equation  at  the  origin.(  See  Eq.(2b).)  Therefore, 
the  before-recruitment  (solid)  curve  lies  above  the  after-recruitment  (dotted)  curve  in  the 
third  quadrant  (p  <  0,  V  <  0)  for  Data  1  for  which  A  (before)  =  0.730  •  A  (after).  For 
Data  6  with  A  (before)  being  close  to  A  (after)  two  cmrves  are  nearly  identical.  In  Data  7 
the  after-recruitment  curve  lies  slightly  above  the  before-recruitment  curve  as  A  (before)  = 
1.143  •  A  (after).  The  data  sets  in  Fig.  11  all  have  the  two  ratios  well  above  unity  with  the 
high  A  ratios  resulting  in  the  after-recruitment  curves  to  lie  above  the  before-recruitment 
curves,  and  the  high  po  ratios  of  Data  13  and  15  being  reflected  in  the  large  shifts  in 
triangles  between  the  two  curves  in  the  p  —  V  diagram. 

From  Eq.(2a)  the  maximum  local  compliance  (=  dV/dp  atp  —  po)  may  be  expressed  as 


dV  ,  ;  A  •  AV 

^(p=po)  = 


(7) 


Fig.  12  is  a  plot  of  [dV/dp(p  =  po)  after  recruitment  maneuver]  vs  [dV/dp(p  =  po)  before 
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recruitment  maneuver]  of  Data  Somce  B.  Changes  in  the  maximum  local  compliance  (= 
compliance  at  the  inflection  point)  are  small  between  before-  and  after-recruitment  data 
with  a  maximmn  change  of  less  than  0.035  [L/cmH20].  It  should  be  mentioned  again  that 
the  local  comphance  hke  other  parameters  may  be  obtained  mathematically  in  a  continuous 
p-V  model  equation,  once  the  parameters  of  the  equation  are  determined  for  a  specifled 
p-V  cmve. 

Depicted  in  Fig.  13  are  p-V  cmves  and  the  corresponding  p  —  V  curves  of  four  rep¬ 
resentative  data  sets  of  patients  with  ARDS,  drawn  over  their  measurement  ranges.  The 
non-dimensional  p  —  V  curves  in  Fig.l3(b),  which,  we  believe,  are  more  useful  for  data 
examinations  and  interpretations,  yield  the  following  observations: 

(1.)  The  magnitude  of  A,  which  is  represented  by  the  slope  of  ap—V  curve,  is  the  largest 
for  Data  20,  and  the  smallest  for  Data  4. 

(2.)  Since  the  origin  of  a  p  —  V  curve  is  the  location  where  the  local  compliance  is  the 
maximum  (i.e.  p  =  po),  Data  4  and  17  extend  well  into  the  region  of  decreasing  compliance, 
while,  the  pressure  range  of  Data  13  is  limited  to  the  region  of  increasing  compliance. 

(3.)  At  p  =  20  cmH20  (shown  as  a  triangle),  the  compliance  is  still  increasing  for  Data  13 
and  20,  close  to  the 'maximum  for  Data  17  and  decreasing  for  Data  4. 

(4.)  The  two  volmne  asymptotes,  Vu  and  Vl,  are  transformed  respectively  into  V  =  1.0 
and  —1.0;  hence,  the  volume  range  of  Data  13  is  closer  to  the Tower  asymptote,  while, 
the  overall  volume  change  of  Data  4  is  small  compared  to  AV  (difference  between  the 
asymptotes). 

Furthermore  the  magnitude  of  A  in  Table  1  indicates  that  there  is  no  lower  corner 
pressme  for  Data  4.  For  each  inflation  data  set.  Table  1  lists  the  maximum  local  compli¬ 
ance.  Its  value  ranges  between  0-03  and  0.11  [L/cmH20]]  much  smaller  values  compared 
to  the  data  from  healthy  adults  shown  in  Fig.l2.  The  maximum  local  compliance,  as 
shown  in  Eq.(7),  is  proportional  to  the  product  of  A  and  AV,  and  inversely  proportional 
to  pressmre  at  the  inflection  point.  Since  both  A  and  po  are  roughly  in  the  same  order  of 
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magnitude  as  healthy  adults  (Fig.8),  the  factor  contributing  most  to  smaller  values  of  the 
Tnavimmn  local  compliance  for  patients  with  ARDS  is  AV  as  evidenced  by  its  values  listed 
in  Table  1  and  2. 

Summary  and  Conclusions 

The  sigmoidal  (tangent  hyperbolic)  equation  is  known  to  represent  various  quasi-static 
p-V  curves  very  closely  [11]-  In  the  present  study  it  has  been  shown  that  the  error  function 
model  equation  also  represents  quasi-static  p-V  curves  well  (Figs.4,5).  Major  parameters 
of  both  the  sigmoidal  (tangent  hyperbolic)  and  the  error  function  model  equations  are 
the  non-dimensional  compliance,  A,  the  maximum  local  compliance,  po,  the  upper  (or 
lower)  volume  asymptote,  Vu  (Vl),  and  the  maximum  voliune  available  for  inflation,  AV. 
Although  both  continuous  model  equations  are  antisymmetric  with  respect  to  po,  the  non- 
dimensional  parameter.  A,  as  well  as  two  volmne  asymptotes  slightly  differ  between  the 
two  equations  as  those  function-specific  parameters  are  selected  to  follow  a  specified  p- 
V  curve  as  closely  as  possible  (Fig.  6).  Two  inflation  data  soiuces,  patients  with  ARDS 
(Data  Source  A)  and  healthy  adults  (Data  Source  B),  are  analyzed  in  detail  using  the  error 
function  p-V  model  equation  with  the  following  results; 

1.  The  alveolar  recruitment  maneuver  lowers  the  pressure  at  the  maximum  compliance, 
Po;  while,  A  remains  roughly  the  same  or  decreases  in  magnitude  (Figs. 8(a),  9).  It  also 
reduces  the  upper  volume  asymptote,  Vu,  substantially  (Fig.8(b)).  The  combined  effects 
of  these  parametric  changes  due  to  the  maneuver  extend  the  range  of  p-V  curves  after  the 
maneuver  further  into  the  region  beyond  the  location  of  the  maximum  compliance  (Figs.  10, 
11). 

2.  The  range  of  po  and  A  covered  by  the  patients  with  ARDS  is  wider  than  the  corre¬ 
sponding  ranges  of  healthy  adults  (Fig.8(a)).  Substantially  lower  magnitudes  of  the  upper 
volume  asymptote,  Vu,  and  the  actual  volume  change  (Fig.8(b))  result  in  lower  values  for 
the  maximum  local  compliance  (Table  1)  compared  to  that  of  healthy  adults  (Fig.  12). 

3.  The  non-dimensional  p  —  V  cxnves  combined  with  the  magnitudes  of  the  four  parame- 
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ters  of  the  model  equation  are  shown  to  help  understand  quantitatively  the  effects  of  the 
recruitment  maneuver  as  well  as  differences  among  patients  with  ARDS  (Figs.ll,  13). 

4.  An  important  advantage  of  the  continuous  model  equations  is  that  various  parameters 
characterizing  the  shape  and  the  range  of  p-V  ciuves,  such  as  the  maximum  local  compli¬ 
ance,  the  pressiue  at  the  maximum  local  compliance  and  the  upper  and  the  lower  corner 
pressures,  may  be  evaluated  readily  from  the  model  equation  once  the  parameters  of  the 
equation  are  determined  from  p-V  curve  data.  The  mathematically  exact  relations  among 
the  parameters  also  implies  that  the  magnitude  of  either  pressure  or  volume  at  a  certain 
location  along  a  p-V  curve  may  only  be  interpreted  correctly  when  compared  to  a  char¬ 
acteristic  pressure  or  volume  of  the  p-V  curve,  as  demonstrated,  for  example,  in  Figs.ll 
and  13  when  we  discussed  the  location  of  p  =  20  cm/faO  relative  to  the  presstue  at  the 
maximum  local  compliance. 
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Captions  for  Tables  and  Figures 

Table  1.  Parameters  of  Error  Function  Equation  for  Data  Source  A 

Table  2.  Parameters  of  Error  Function  Equation  for  Data  Source  B 

Fig.  1.  Continuous  p-V  model  equation. 

Fig.  2.  Non-dimensional  pV  curve  corresponding  to  Fig.l. 

Fig.  3.  Example  of  p-V  curve  from  Data  Source  A.  filled  circle  =  Data  D, 
solid  =  error  function  p-V  equation,  dotted  =  sigmoidal  p-V  equation. 

Fig.  4.  Example  of  p-V  curve  from  Data  Source  B.  Dotted  ciirve  =  piecewise  continuous 

equation,  Eq. (6),  Solid  ciuve  =  error  function  equation,  Eq. (2a), 

Triangle  =  upper  and  lower  data  limits. 

Fig.  5.  (V  -  Vl)/AV  vs  (v/i/4)A{p/po  -  1)  of  Data  Source  a.  .  .  . 

Unfilled  circle  =  inflation.  Filled  circle  =  deflation. 

(  ■  , 

Fig.  6.  (a)  As  (sigmoidal  model  equation)  vs  Ae  (error  function  model  equation), 

(b)  pos  (sigmoidal  equation)  vs  poe (error  function  equation), 
for  inflation  data  sets  of  Data  Source  A. 

Fig.  7.  (a.)  plip  [cmHaO]  vs  Pmci  [cmH20] 

(Allied  (unfilled)  circle  =  before  (after)  recruitment  maneuver). 

(b.)  puip  [cmH20]  vs  pw  [cmH20]  of  Data  Source  B. 

(unfilled  (filled)  circle  =  before  (after)  recruitment  maneuver). 

(p^ct  and  Pmcd  evaluated  from  error-function  equation.) 

Fig.  8.  (a.)  A  vs  po,  (b.)[Vy  -V(p  =  0)]  vs  [V(p  =  20cmIl2O)  -V(p  =  0)] 
for  inflation  data  sets  from  Data  Source  A  and  B. 

Square  =  patients  with  ARDS,  Cross  =  healthy  adults  before  recruitment  maneuver, 
Triangle  =  healthy  adults  after  recruitment  maneuver. 

Fig.  9.  Po  (before  maneuver) /pq  (after  maneuver))  vs  A  (before) /A  (after) 

Two  numbers  are  Data  Set  No.,  followed  by  his  or  her  age  in  the  bracket. 
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unfilled  circle  =  male,  filled  circle  =  female. 

Fig.  10.  p  —  V  curve  and  the  corresponding  p-V  curve  of  error  function  model  equation. 
Data  No.l  (top),  No.6  and  No.7  (bottom)  of  Data  Source  B  (healthy  adults). 

Solid  =  Before  recruitment  maneuver.  Dotted  =  After  recruitment  maneuver. 
Triangle  =  location  of  p  =  20  [cmH20]. 

Fig.ll.  p—V  curve  and  the  corresponding  p-V  curve  of  error  function  model  equation. 
Data  No.  11  (top),  13  and  15  (bottom)  of  Data  Source  B. 

Solid  =  Before  recruitnient  maneuver.  Dotted  =  After  recruitment  maneuver. 
Triangle  =  location  of  p  =  20  [cmH20]. 

Fig.l2.  dV/dp{p  =  Po)  (maximum  local  compliance)  before  recruitment  maneuver 
vs  dV/dp  (p  =  Po)  after  recruitment  maneuver  of  Data  Source  B. 

Fig.l3.  (a)  p  —  V  curve  and  (b)  the  corresponding  pV  curve  of  error  function  model 
equation  for  Data  No.4,  13,  17  and  20  of  Data  Source  A  (patients  with  ARDS). 
Triangle  =  location  of  p  —  20  [cTnH20]  • 
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Table  1.  Inflation  Peirameters  of  Error  Function  Equation  for  Data  Source  A 


Data  No 

A 

Po 

Ay 

(dV/dp)atpo 

1. 

2.9652, 

22.39, 

2.3726, 

-0.0709 

0.0785 

2. 

2.7173, 

21.98, 

1.5612, 

-0.0762 

0.0482 

3. 

3.3532, 

25.08, 

1.6193, 

-0.0365 

0.0541 

4. 

1.6273, 

13.32, 

3.1567, 

-0.5727 

0.0963 

5. 

2.7160, 

30.36, 

1.6216, 

-0.0768 

0.0362 

6. 

2.6029, 

23.88, 

1.5066, 

-0.0962 

0.0410 

7. 

1.7288, 

20.15, 

3.0989, 

-0.3645 

0.0664 

8. 

1.8901, 

14.95, 

1.7905, 

-0.2185 

0.0565 

9. 

3.5379, 

25.24, 

2.7570, 

-0.0350 

0.0965 

10. 

2.7981, 

26.20, 

3.7326, 

-0.1583 

0.0996 

11. 

2.4296, 

17.89, 

1.3424, 

-0.0997 

0.0455 

12. 

1.2500, 

11.59, 

1.2470, 

-0.2829 

0.0336 

13. 

3.2725, 

29.86, 

3.9075, 

-0.0861 

0.1070 

14. 

2.7915, 

15.31, 

1.6284, 

-0.0705 

0.0742 

15. 

1.9487, 

18.37, 

1.8797, 

-0.2423 

0.0498 

16. 

2.4566, 

26.98, 

1.3277, 

-0.0787 

0.0302 

17. 

2.1027, 

19.31, 

1.3316, 

-0.1279 

0.0362 

18. 

1.1672, 

13.98, 

2.0685, 

-0.5149 

0.0431 

19. 

3.1381, 

26.80, 

3.1306, 

-0.0851 

0.0916 

20. 

5.4709, 

30.03, 

1.7695, 

-0.0097 

0.0805 

21. 

3.2818, 

24.43, 

2.8956, 

-0.2075 

0.0972 

Po  in  [cmH20\,  AV  and  Vl  in  [L],  dV/dp  in  [L/cmH20]. 
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Table  2.  Parameters  of  Error  Function  Equation  for  Data  Source  B 
Data  No  A  Po  ^ 


l.F 

B. 

1.5136 

A. 

2.0711, 

2.F 

B. 

2.5520 

A. 

1.7587 

3.F 

B. 

2.0979 

A. 

2.2960 

4.M 

B. 

2.1491 

A. 

2.1887 

5.F 

B. 

2.2192 

A. 

1.5729 

6.M 

B. 

2.3950 

A. 

2.4387 

7.M 

B. 

2.6997 

A. 

2.3615 

8.F 

a 

2.2901 

A. 

1.9974 

9.F 

B. 

— 

A. 

2.0340 

lO.M 

B. 

3.4081 

A. 

2.3998 

ll.M 

B. 

3.1746 

A. 

2.2828 

12.F 

B. 

2.5060 

A. 

2.3844 

13.F 

B. 

2.3600 

A. 

1.8107 

14.M 

B. 

2.2260 

A. 

2.0800 

15.F 

B. 

3.5047 

A. 

2.1561 

25.00 

14.48 

5.1043 

2.7130 

31.70 

18.17 

5.4303 

4.1416 

21.76 

18.19 

3.7676 

3.1331 

21.11 

18.57 

4.5443 

4.4023 

24.60 

13.46 

6.4953 

5.2884 

23.21 

19.40 

4.2092 

3.9049 

31.00 

16.85 

7.3180 

5.5288 

30.80 

17.20 

8.2489 

5.0638 

15.90 

3.7392 

28.80 

17.92 

7.211- 

5.7276 

25.86 

19.70 

7.0726 

6.79820 

23.20 

16.36 

5.0810 

4.3799 

29.17 

14.00 

6.81755 

4.03567 

24.55 

18.22 

6.7170 

5.5875 

25.70 

14.65 

4.5022 

3.3016 

-0.8889 

-0.2700 

4.2154 

2.4430 

-0.2863  . 

-0.6590 

5.7166 

3.4826 

-0.3399 

-0.2300 

3.4277 

2.9031 

-0.3726 

-0.3554 

4.1717 

4.0469 

-0.4971 

-0.8895 

5.9982 

4.3989 

-0.2567 

-0.2722 

3.9525 

3.6327 

-0.2167 

-0.5641 

7.1013 

4.9647 

-0.5185  . 

-0.5478 

7.7304 

4.5160 

-0.4068 

3.3324 

-0.0491 

-0.5502 

7.1627 

5.1774 

-0.1361 

-0.8620 

6.9365 

5.9362 

-0.2476 

-0.3.799 

4.8334 

4.0000 

-0.4353 

-0.6883 

6.3822 

3.3472 

-0.4828 

-0.6761 

6.2342 

4.9114 

-0.0027 

-0.4836 

4.4995 

2.8180 

B  =  Before  recruitment  maneuver,  A  =  After  recruitment  maneuver. 
Po  in  [cmH20],  AV,  Vi  and  Vy  in  [L]. 
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Fig.1 .  Continuous  p-V  model  equation. 
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Fig.2.  Non-dimensional  p-V  curve  corresponding  to  Fig.  1. 
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Fig. 3.  Example  of  P-V  curve  from  Data  source  A. 
filled  circle  =  Data  D,  solid  =  error  function  p-V  equation,  dotted  =  sigmoidal  p-V 
equation 


Fig.4.  Example  of  p-V  curve  from  Data  Source  B. 
tted  curve  =  piecewise  continuous  equation,  Eq(6),  Solid  curve  =  error  function  equation, 
r2a).  Triangle  =  upper  and  lower  data  limits. 


Fig.5.  (y-V^)I^V  vs  (^/^ IA)K{pl po  -1)  of  Data  Source  A. 
Unfilled  circle  =  inflation,  Filled  circle  =  deflation. 


0  id  20 

Po,s 

Fig.6.  (a)  Ag  (sigmoidal  model  equation)  vs  A,  (error  function) 
model  equation.(b)  ^  (sigmoidal  equation)  vs  ^  (error  function 
equation)  for  inflation  data  sets  of  Data  Source  A. 


30 


30 


Pmcd 


Fig.7.  (a)  [cmH20]  vs  p^,.[cmH20](filled  (unfilled)  circle=  before 

(after)  recruitment  maneuver),  (b)  /7y,p[cmH20]  vs  [cmH20]  of  Data 
Source  B. (unfilled  (filled)  circle  =  before  (after)  recruitment  maneuver, 
and  p„^^  evaluated  from  errorfunction  equation.) 
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V(P=20)-V(P=o)  [L] 

Fig.8.  (a.)  A  vs  p,  ,  (b.)  [V^  -V(p  =  0)]  vs  [Vip  =  20cmH,0)-V{p  =  0)] 
for  inflation  data  sets  from  Data  Source  A  and  B. 

Square  =  patients  with  ARDS,  Cross  =  healthy  adults  before  recruitment  maneuver, 
Triangle  =  healthy  adults  after  recruitment  maneuver. 
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Fig.  9.  Po  (before  maneuveii after  maneuver  vs  A  (before)/  A  (after). 

Two  numbers  are  Data  Set  No.,  followed  by  his  or  her  age  in  the  bracket, 
unfilled  circle  =  male,  filled  circle  =  female. 


V  ' 

,  0.5 

.A'* 

5  -1  ( 

■IHH 

-0.5 

-1 

■■ 

Data  No.  1 


Data  No.  6 


P[cmH20| 


Data  No.  7 

Fig.10.  p-V  curve  and  the  corresponding  p-V  curve  of  error  function  model  equation. 
Solid  ^Before  recruitment  maneuver,  Dotted  =  After  recruitment  maneuver,  Triangle  =  location 
of  p=20  [cmH20]. 
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=  Po)  (maximum  local  compliance)  before  recruitment  maneuver 
recruitment  maneuver  of  Data  Source  B. 
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Chapter  2.  A  Mechanistic  Model  for  Quasi-Static  Pulmonary  Pressure- Volume 
Cmves:  Model  Development  for  Inflation  Process 

Abstract 

A  mechanistic  model  of  a  total  respiratory  system  is  proposed  to  understand  differences 
in  quasi-static  pressure-volume  (p  —  V)  curves  of  the  inflation  process  in  terms  of  the  alve¬ 
olar  recruitment  and  the  elastic  distension  of  the  wall  tissues.  In  the  model,  based  on  the 
Boltzmann  statistics,  the  total  respiratory  system  consists  of  a  large  number  of  elements, 
each  of  which  is  a  subsystem  of  a  cylindrical  chamber  fitted  with  a  piston  attached  to 
a  spring.  The  alveolar  recruitment  is  simulated  by  allowing  a  distribution  of  the  critical 
pressure  at  which  an  element  opens;  while  the  wall  distension  is  represented  by  the  piston 
displacement.  Various  parameters  in  the  error-function  p  —  V  model  equation  are  related 
to  the  properties  of  the  mechanistic  model.  The  parameters  of  the  model-based  p-V  equa¬ 
tion  are  determined  for  each  clinical  data  set  for  a  total  of  twenty  one  p  —  V  data  sets  of 
patients  with  acute  respiratory  distress  syndrome  by  a  computational  minimization  proce¬ 
dure  between  the  equation  and  the  data  points,  results  of  which  show  excellent  agreements 
between  the  two. 
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Nomenclature 


A 

As 

B 

C 

f,  F 

k 

N 

Nj 

P 

P 

Pc  3 
Pc  3 

Pf 

Po 

Uip) 

AU 

V 
Vp 

Vl(u) 

V 

% 

AV 

Vb 

Vj 


{k/As)yo  =  Po 

piston  surface  area  on  which  pressure  is  acting  (Fig.4) 

{k/Aa)yT  =  Po  •  yro)  threshold  pressure  for  onset  of  saturation 

(^/2)i/2.po  =  \/^A/4 

distribution  functions  (Eqs.(6)  &  (13c)) 

functions  defined  in  Eq.(13b) 

spring  constant  [N/m]  (Fig.4) 

total  number  of  TRS  elements 

number  of  elements  at  energy  level  j 

pressmre  (interpleiural  pressure  difference) 

non-dimensional  pressure,  p/po  —  1 

critical  pressure  at  which  an  element,  j,  ‘pops  open’. 

Pcj/Po  (Eq.(13c)) 

pressure  at  the  end  of  inflation  process 

pressTure  at  the  inflection  point  in  p-V  equation,  Eq.(l) 

total  energy  of  TRS  at  p=p 

=  U(p  =  p,)  -  U(p  =  0) 

volume 

volume  change  from  the  state  of  p  =  0 
lower  (upper)  volume  asymptote  (Fig.l) 

non-dimensional  volume,  (V  —  ((Vu  +  Vl)/2)/(AV/2),  (Eq.(2)) 
volume  of  an  element  j 
Vu  —  Vl  =  NVoiyro  +  1) 

‘pop-open’  volume  (=  AsPo)  (Fig.4) 
piston  displacement  of  an  element  j  (Fig.4) 
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.yN,. 

yo  ‘pop-open’  displacement,  (=  Vq/As)  (Fig.4) 

yx  piston  stroke  limit  (Fig.4) 

Vto  Vt/vo 

Greek  symbols: 

constant  of  proportionality  (Eq.(la)) 
constant  in  the  distribution  function  (Eq.(4)) 
=  (Al/k)  (3  {Eq.{5)) 
energy  stored  in  an  element  j  (Eq.(3)) 
opo A V (non-dimensional  parameter)  (Eq.(l)) 
(8/7r)^/^/A,  Standard  deviation  (Eq.(13c)) 
Ap/2  (Eq.(2)) 


a 

p 

p 


^0 


u> 


Superscript: 


^  related  to  a  single  TRS  element 

Acronyms: 

ARDS  acute  respiratory  distress  syndrome 

TRS  total  respiratory  system 
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Introduction 


Quasi-static  pulmonary  p  —  V  (pressure  -  volume)  curves  are  used  routinely  to  obtain 
quantitative  information  on  the  respiratory  system  that  is  important  for  both  research  and 
clinical  guidances,  as  the  conditions  of  gas  exchange,  the  primary  role  of  the  respiratory 
system,  are  related  to  the  characteristics  of  the  cmve.  Dming  the  inflation  (inspiration)  and 
the  deflation  (exhalation)  processes,  the  respiratory  system  changes  its  volume  (measured 
in  L  (=  10“^  m^)  or  mL),  lung  (alveolar)  pressure  as  well  as  the  pleural  pressure  (the 
pressure  of  the  thin  liquid  film  that  couples  the  lungs  and  the  chest  wall  pleurae).  The 
pressure,  p,  refers  to  the  interpleural  pressme  difference  (i.e.  the  difference  between  the  lung 
pressure  and  the  pleural  pressure)  measured  in  water  head  [cm-H20]  (1  cmH20  =  98  Pa). 
Clinical  p-V  curves  are  commonly  obtained  for  an  anesthetized  human  subject  in  supine 
position  by  sequentially  adding  (or  withdrawing)  incremental  gas  volmnes  (~  50-100  mL) 
in  a  stepwise  manner  (with  a  duration  of  ~  5  seconds  per  step)  [1,2].  Fig.l  is  a  t3q)ical 
inflation  p-V  cmve,  consisting  of  a  nearly  linear  region  of  high  compliance  (i.e.  large 
dV/ dp)  sandwiched  between  two  segments  with  low  compliance  at  low  and  high  pressure 
regions.  The  shape  of  the  p  —  V  cmve  is  affected  by  two  mechanisms,  the  distension  of  the 
elastic  respiratory  wall  tissue  components  and  the  recruitment  (‘pop-open’  mechanism) 
of  the  alveoli.  The  latter  is  the  opening  of  alveoli  overcoming  the  surface  tension  at 
the  interface  between  the  gas  and  the  liquid  film  lining  the  alveolar  surface.  A  pressure 
increase  (i.e.  an  increase  in  the  interpleural  pressure  difference)  results  in  the  recruitment 
of  a  greater  number  of  alveoli.  The  high  compliance  is  believed  to  be  associated  with  both 
the  distension  of  open  parts  of  the  lungs  and  the  (alveolar)  recruitment  of  collapsed  parts 
of  the  lungs  [3].  Some  protective  ventilation  strategies,  based  on  patients’  quasi-static  p-V 
cmves,  have  been  proposed  for  Irmg  disease  patients  in  intensive  care  miits.  Amato  and 
coworkers  [4,5]  demonstrated,  based  on  their  clinical  study  involving  patients  with  acute 
respiratory  distress  syndrome  (ARDS),  that  a  ventilator  strategy  guided  by  the  p)-V  curve 
resulted  in  reduced  lung  trauma,  a  high  weaning  rate  and  improved  smvival  compared  with 
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a  conventional  ventilator  strategy  without  the  p-V  curve  guidance.  Also,  a  recent  ARDS 
Network  report  [6]  on  a  clinical  study  involving  861  patients  shows  lower  mortality  in  the 
group  treated  with  lower  tidal  volume  than  in  the  group  treated  with  traditional  higher 
tidal  volumes.  Although  a  use  of  p-V  curves  is  not  mentioned,  the  report  underscores  the 
importance  of  optimized  ventilator  strategy. 

In  order  to  quantify  the  characteristics  of  p-V  curves  as  well  as  their  changes  observed 
in  clinical  settings,  various  p-V  model  equations  have  been  proposed  [7-12].  Parameters  in 
model  equations  are  determined  from  statistical  processing  of  chnical  data.  It  is  important 
that  these  parameters  should  have  some  physiological  interpretations.  Also,  to  understand 
the  shape  of  p-V  curves  in  terms  of  mechanical  behavior  of  lungs,  multi-compartment  lung 
models  were  developed  and  used  to  obtain  information  on  the  effects  of  lung  elasticity  and 
a  degree  of  alveolar  recruitment  on  p-V  ciuves  [3,13].  Although  these  analyses  serve  to 
relate  the  internal  elastic  conditions  of  the  total  respiratory  system  (TRS)  to  general  p-V 
cmve  behavior,  there  has  been  no  attempt  to  interpret  individual  differences  in  p-V  curves 
directly  in  terms  of  internal  elastic  properties,  alveolar  recruitment  and  their  changes. 
Prom  an  analytical  viewpoint  the  quasi-static  p-V  ciuves  are  more  amenable  to  theoretical 
studies  because  at  each  state  we  may  be  able  to  apply  equilibrium  principles.  An  overall 
objective  of  this  report  is  to  test  the  hypothesis  that  a  mechanistic  model,  based  on  the 
continuous  alveolar  recruitment  and  the  elastic  distension  of  the  wall  tissues,  is  effective  in 
understanding  a  relation  between  the  observed  pulmonary  behavior  (as  p-V  curves)  and 
the  corresponding  internal  respiratory  response  (in  terms  of  the  mechanistic  model). 
Continuous  Equation  for  Quasi-Static  p-V  Curves 

In  the  past  piecewise-continuous  equations  were  used  to  generate  such  quantities  as 
comphance,  the  lower  and  upper  inflection  points  that  may  reflect  the  internal  conditions 
of  TRS  [7-12].  There  are  two  continuous  model  equations  that  simulate  various  p)-V  curves 
accurately  over  the  entire  range  of  p-V  data.  One  is  a  sigmoidal  (tangent  hyperbolic) 
equation,  and  the  other  an  error-function  equation,  both  originally  proposed  by  Venegas 
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and  his  coworkers  [1].  Since  the  analytical  development  to  follow  utilizes  the  error  function 
representation  of  p-V  curves,  its  characteristics  are  discussed  below  in  some  detail. 

The  error  function  p-V  equation,  plotted  in  Fig.l,  may  be  expressed  as 


where  AV  =  Vu  -  Vl,  Vu  =  the  upper  asymptote,  Vl  -  the  lower  asymptote,  a  = 
positive  constant,  po  =  &  pressure  at  the  midpoint  (inflection  point)  of  the  curve  and  A 
(non-dimensional)  =  apoAV  [1,2,14]  with 

erf{x)  =  -^f  e~^^dt  er/(oo)  =  1,  erf{-x)  =  -erf{x). 

Jo 

The  corresponding  non-dimensional  forms  are, 

^  =  f  V  =  erf{^i^).  (2) 

where  V  =  [V  —  {Vu  +  Vl)/2]/(AF/2),  p  =  (p/po)  —  1,  oj  =  Apl2. 

Eq.(2)  satisfies  the  following  conditions: 

V  {p  =  0)  =0,  y  ip)  =  -y  (—p)  (antisymmetry  with  respect  to  p  =  po)) 

dF/dp  (p  =  0)  =  A/2,  F  (p  ±oo)  =  ±  1,  dF  (p  ±oo)/dp  =  0. 

A  clinical  data  somrce  of  p-V  cmrves  we  use  in  the  present  analyses  are  twenty  one  data 
sets  of  ARDS  patients  (both  inflation  and  deflation  data)  in  supine  position  by  Harris  et 
al  [2],  made  available  by  the  authors.  Data  sets  are  analyzed  by  minimizing  the  difference 
between  data  points  and  the  error  function  model  equation  through  the  application  of  the 
method  of  least  squares  to  obtain  the  parameters,  A,AV,po  and  Vy  (or  Vi,).  Plotted  in 
Fig.2  (Fig.3)  are  264  inflation  data  points  (225  deflation  data  points)  in  terms  of  Eq.(lb), 
(V  -  Vl)/AV  vs  0rA(p/po  —  l)/4.  Agreement  is  excellent  between  the  data  and  the  p-V 
equation  with  B?  (the  coefficient  of  determination)  =  0.99938  for  the  inflation-  and  = 
0.99907  for  the  deflation-  data  points. 
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Development  of  a  Mechanistic  Model 

An  overall  objective  of  the  development  of  a  mechanistic  TRS  model  is  to  derive  an 
equation  for  the  volume  participating  in  the  p-V  variations.  We  consider  a  TRS  comprised 
of  a  very  large  nmnber  of  elements  with  N  =  total  number  of  elements.  Based  on  the  char¬ 
acteristics  of  a  single  element  that  are  common  to  all  elements,  a  distribution  of  elements 
is  derived  over  a  distribution  parameter.  The  mechanistic  model  of  an  element  is  shown 
in  Fig.4.  An  arbitrary  element,  j,  consists  of  a  cylindrical  chamber  containing  a  piston 
(with  its  surface  area,  Ag,  which  is  attached  to  a  spring  with  its  spring  constant, 

k  [N/m].  The  element  is  closed  when  the  piston  is  located  at  the  left  end  of  the  cylinder  in 
Fig.4.  When  pressure  acting  on  the  left  end  of  the  piston  reaches  a  certain  critical  value, 
the  piston  suddenly  moves  to  a  new  position  (‘pop-open’  mechanism)  with  the  elemental 
volume,  V^,  in  the  figure  indicating  an  elemental  volume  increase  due  to  the  sudden  piston 
displacement  of  yo-  Once  the  element  is  open  with  its  volume  of  Vo?  any  further  increase  in 
pressure  results  in  a  volume  increase  as  the  piston  moves  to  the  right  until  it  reaches  the 
end  of  the  cylinder.  (The  symbol,^,  indicates  an  elemental  quantity.)  In  the  model  the 
pop-open  volume,  Vq  {—  Ag  yo),  and  a  further  volume  increase  due  to  piston  displacement 
represent  the  opening  of  alveoli  and  the  elastic  distension  of  the  wall  tissues  respectively. 
We  define  Pcj  as  the  critical  pressure  at  which  the  element,  j,  ‘pops  open’. 

y«s. 

Referring  to  Fig.4,  the  elemental  volume,  Vj,  at  p  (>  Pcj)  is  equal  to  Vo  +  AaPj',  which, 
upon  application  of  a  quasi-static  force  balance  across  the  piston,  Ag  (p—Pcj)  =  kyj,  may 
be  expressed  as,  •+■  (A^/fc)  {p-Pcj)-  Also  the  piston  position  of  an  element  reaches 

its  stroke  limit  of  yr  when  pressure,  p,  reaches  {j>cj  +  {k/Ag)yT)-  The  mechanistic  model  of 
an  element,  therefore,  goes  through  three  stages  in  the  inflation  process  —  closed,  open  & 
unsatmated  (i.e.  yj  <  yr)  and  open  &;  saturated  (i.e.  yj  =  yr)-  The  model  assumes  that 
the  critical  ‘pop-open’  pressure,  pcj,  as  well  as  the  location  of  the  piston  for  open  elements, 
yj,  vary  fi:om  element  to  element  at  an  arbitrary  quasi-static  state  {p,V),  and  that  other 
quantities  such  as  k,  Ag,  %  and  yr  are  constant  and  common  for  all  elements.  The  energy 
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level  of  an  open  and  unsatiu-ated  element,  j,  consists  of  the  activation  energy  required  to 
pop  open  the  element,  tjA,  and  the  energy  stored  in  the  spring,  ejs^  For  ZjA,  we  assign 
the  compression/expansion  work  under  constant  pressure;  i.e.  tjA  =  PcjVo'i  while,  ejs  is 
equal  to  k  ^|/2,  which  may  be  expressed  in  terms  of  pc  j  astjs  =  {Al/2k){p  -  from 
an  application  of  the  force  balance. 

In  summary,  at  a  quasi-static  (-equilibrium)  state  at  p  =  p,  a  TRS  element,  j,  belongs  to 


one  of  the  following  states: 

if  p  <  Pcj,  the  element,  j,  is  closed  with  Vj  =  0,  tj  =  0. 

if  p^j  =  p,  the  element,  j,  pops  open  with  Vj  =  Vo,  tj  =  Pc jVo. 

if  p  -  {k/As)yT  <  Pcj  <  P,  the  piston  of  the  open  and  unsaturated  element,  j, 
moves  to  a  location,  pj,  with  Vj  =  Vq  +  Ag  Vj,  tj  =  PcjVo  +  (At/2k){p  —  Pcj)^- 
if  p^j  <p  —  {klAs)yT,  the  piston  of  the  open  and  saturated  element,  j, 

remains  at  the  stroke  limit,  yT,  with  Vj  =  Vq  +  Ag  Pr,  —  PcjVo  +  [k/2)y^. 

The  state  of  an  element  follows  the  sequence  above  during  the  inflation  process  as  p  in¬ 
creases.  To  obtain  an  explicit  form  of  the  distribution  function  of  TRS  elements  over  the 
distribution  parameter,  Pcj,  we  focus  on  open  and  unsaturated  elements,  for  which  the 
elemental  energy  shown  above  may  be  rewritten  as. 


tj  {open,  unsaturated)  = 


(3) 


According  to  the  Boltzmann  statistical  model  ([15,16]  for  example),  which  assumes  that 
there  is  no  limit  in  the  number  of  elements  per  energy  state,  the  most  probable  distribution 
Nj/N  (a  fraction  of  elements  at  an  energy  level,  tj),  may  be  expressed  as 


Nj  /  N  =  —  unspecified  constant)  (4) 

3 

A  substitution  of  Eq.(3)  into  Eq.(4)  with  the  summation  replaced  by  an  integral  over  the 
whole  range  of  Pcj  for  a  large  number  of  elements,  yields 
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where  =  (A^k)  0.  It  should  be  noted  that  the  integration  in  the  denominator  ranges 
from  -  oo  and  +  oo.  As  summarized  above,  any  arbitrary  element  remains  active  (open 
and  unsaturated)  only  in  a  certain  range  of  Pcj]  therefore,  the  application  of  Eq.(5)  over 
the  entire  range  of  —  oo  <  Pc  j  <  oo  assumes  that  the  distribution  function  that  is  valid  for 
active  elements  is  also  applicable  in  evaluating  the  number  of  closed  as  well  as  saturated 
elements.  Then,  upon  performing  the  integration  of  the  denominator  in  Eq.(5),  we  obtain 
for  dNj/N{—  a  number  fraction  of  elements,  for  which  the  magnitude  of  pc  j  ranges  between 
Pcj  a.nd  Pc j +dpcj), 


Noting  that  elements,  j,  with  pcj  in  the  range  of  0  <  Pcj  <  p,  are  open  at  p  =  p,  and 
that  (k/As)  Vt  (=  B)  is  the  pressure  at  which  an  element  j  with  Pcj  =0  reaches  the  piston 
stroke  limit  of  yr,  the  volume  change  with  pressure  needs  to  be  evaluated  for  the  following 
two  pressure  ranges;  pressure  range  1:  0  <  p  <  JB  and  pressure  range  2:  B  <p. 

Pressure  Range  1:  0  <  p  <  B 

Since  the  pressmre  is  below  B  (the  threshold  pressure  for  the  onset  of  saturation),  all 
open  elements  are  active  (unsaturated)  with  pj  <  yr-  Then  the  total  volume,  Vp{=  V{p  = 
p)-V{p  =  0)-,  i.e.  a  volume  change  from  the  state  of  p  =  0),  participating  in  the  inflation 


process  is. 


»p  ^  poo 

Vp{p)  =  N  /  {Vo  +  Asyj)f{p  =  p)dpcj  +  Asyj  /  f{p  =  p)dpcj  . 
[7pc,=o  Jpc}=p 

The  first  term  on  the  right  hand  side  represents  a  volume  increase  due  to  the  elements  that 
pop-open,  followed  by  a  piston  displacement  (yj  —  A\{p  —  Pcj)/k  from  the  force  balance) 
during  the  inflation  process  from  p  =  0  to  p  =  p ;  while,  the  second  term,  noting 

poo  poo 

/  f{p  =  p)  dpcj  =  f{p  =  0)  dpcj , 

Jpcj=P  •'Pcj=0 

accounts  for  the  elements  that  are  already  open  at  p  =0  and  the  piston  displacement 
yj  (=  A}gp/k)  is  the  only  mechanism  available  for  the  volume  increase. 


47 


After  expressing  yj  in  terms  of  p  and  pc  j  as  shown  above,  the  equation  may  be  written  as, 


/p-A  a2  /?  roo  a2 

{2Vo  + -^t)-exp{-^t^)dt  +  -^p  •  exp{- -z^)dz 


(7) 


where  A  —  {k/Ag)  yo,  t  =  —z  =  (p  —  Pcj  —  k  Vo/Ag). 

Pressure  Range  2:  B  <p 

In  this  pressure  range,  the  elements  with  0  <pcj  <  p  —  B  as  well  as  the  elements  that 
are  already  open  at  p  =0  are  saturated  (i.e.  yj  =  yx  for  the  elements);  while,  the  elements 
with  p  —  B  <  Pcj  <  p  remain  unsaturated  (i.e.  yj  <  yr)]  therefore, 

rv-B 


fP-B  ^  '■- 

Vp  =N  I  (Vo  +  Asyx) '  f(p  —  p) dpcj  +  N  I  (Vo  +  Agyj )’  f(p  —  p) dpcj 
Jo  Jp—B 

i*00 

+  N  AsPt  ■  fip  =  p)  dpcj. 

Jp 

Eqs.(7)(8),  after  integration,  become, 

(0  <  p  <  B)  =  wfi,  [l,  +  i  |(1  -  /i)  -  . 

Vp(  B<p  )  =  NVo\h  +  ^(l-h)  +  —[^ 


(8) 


(9o) 


h- 


yro  + 1 


2y/nC 


/sb)]  (9i>) 


where 


li=erf(0,  l2(i>)  =  exp{-CH^-lY)-exp(-C^),  h(p)  =  eTf(C  1)), 
I4  =  er/(C'(l  —  yro)))  h  =  &xp{—C^{yTo  —  1)^)  —  exp{—  C^), 

C=  (;0/2)'^".A,  yTo  =  yT/m  =  B/A. 

Therefore,  the  mechanistic  model  yields  the  following  p-V  equation; 


V  {0  <  p  <  B)  =  Vu 

■il-Ii)jP 
A 


ArO  ( P  ^  \  yTO  d”  f  yro  1  |-  ,  -^5  ^2\P)  ,  f  ( 

+  «Vo[^-^(j-VTo) - 2 - 2 — + 


h  -  h{p) 


.  ^  iVyo(yro  +  l)  NVo{yTo  +  '^)  j  ,  ^ 

V{B  <p  )  =  Vu - - - +  - - ■'alPj- 


(106) 
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In  order  to  relate  the  present  TRS  model,  intended  to  describe  the  internal  respira¬ 
tory  conditions,  to  the  p-V  curve  that  quantitatively  describes  overall  variations  in  TRS 
conditions,  an  additional  condition  that  needs  to  be  satisfied  is  conservation  of  energy.  For 
a  quasi-static  process  fi’om  an  initial  state  of  (p  =  0,  =  V{p  =  0))  to  a  final  state  of 

(p  =  py,  V  =  V(p  =  pf))  {pf  =  the  pressure  at  the  end  of  a  measiured  i>-V  curve),  the 

conservation  of  energy  neglecting  any  dissipative  mechanisms  may  be  written  as, 

-V(p=p/) 

AU{=U{p  =  Pf)  -  U{p  =  0))  =  pdV,  (11) 

Jv(p=o) 

where  U{p)  represents  the  total  energy  of  TRS  at  p  =  p  to  be  evaluated  from  our  mecha¬ 
nistic  model;  while  the  right  hand  side  of  the  equation  is  work  associated  with  the  inflation 
process  that  must  be  evaluated  firom  the  p-V  model  equation.  A  further  development  of 
the  energy  equation  will  be  discussed  in  the  next  section. 

Mechanistic  Model  vs  Error  Function  p-V  Equation 

Relations  between  the  parameters  in  the  error  function  p-V  equation  and  the  pa¬ 
rameters  in  om:  mechanistic  model  are  derived  based  on  the  observation  that  the  p  —  V 
relation,  Eq.(lOb),  for  the  high  pressvure  region  as  well  as  the  corresponding  equation  for 
the  local  compliance. 


V{B<p)  =  Vu- 
{B<  p) 


NVo(yTo  +  1)  ^ 

2  2  A 


,  NVoivTO  +  l)  2C 

^(B<p)  =  — 2 - I’ 


become  identical  to  the  error  function  model  equation,  Eq.(l), 

if  we  set 


NVoivTo  +  1)  =  AV,  a{=  {k/As)yo)  =Po,  C  (=  •  a)  =  ^  A.  (12) 
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Before  further  developments  are  made  on  the  mechanistic  model,  our  results  are  summa¬ 
rized  below,  based  on  the  parametric  relations,  Eq.(12),  between  the  error  function  p-V 
equation  and  the  model. 

Pressure  -  Volume  {p  —  Vp)  Equation: 


Vp{0  <p<po‘  Vto)  = 
Vp  {po  •  Vto  <p)  = 


Vto  + 1 

AV 


Vto  + 1 


+  (13a) 


Pressme  -  Volume  {p  —  V)  Equation: 


V  ( 0  <  p  <  Po  •  yTO )  = 
AV 


Vu  +  Vl 


Vto  + 1 


ttA 


r 


xr/  -  ^  \  Vb-t-Vx,  AV  .  . 

V(po-yTo<P)  = - 2 - 

where  h  =  erf{C),  hip)  =  expi-C^p^)  -  exp(-C^), 
hip)  =  erfiCp),  h  =  er/(C'(l  -  ^)), 
h  =  expi-C"^  ivTo  -  1)^)  -  exp(-  C^),  C  =  V^A/4. 


Distribution  Function: 


(136) 


dNj 

N  •  dpc  j 


where 


Pcj  ^Pcj/POy  ^  (®/^)  ^ • 


(13c) 


The  model-based  p-V  equation,  Eq.(13b),  consists  of  two  regions.  The  solution  for  the 
high  pressure  region  is  identified  with  the  error-function  p-V  equation.  The  p-V  equation 
for  the  lower  pressure  region  as  well  as  the  boundary  pressure  between  the  high-  and 
low-pressure  region  contain  the  parameters  of  the  p-V  equation.  A,  po,  AV,  Vu  (orVc)), 
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and  an  additional  parameter,  yxo-  Conservation  of  energy,  Eq.(ll),  is  utilized  to  find  the 
magnitude  of  yro-  Similar  to  the  p-V  model  equation,  Eq.(13b),  the  evaluation  of  Eq.(ll) 
depends  on  the  magnitude  of  the  final  pressure,  pf,  relative  to  the  boundary  pressure, 
PoVto,  between  the  high  and  the  low  pressure  regions  of  the  model-based  p-V  equation. 
As  will  be  shown  later  in  the  analyses  of  clinical  data,  the  magnitude  of  yTo  is  less  than 
unity;  hence,  the  conservation  of  energy  is  further  developed  for  the  case  of  po‘2/ro  ^  Pf  (i-G- 
yro  —  1  <  Pf{—  PfIPo  —  !))•  Accordingly  the  left  hand  side  of  Eq.(ll)  may  be  evaluated 
from  the  elemental  distribution  function,  F{p)  of  Eq.(13c),  along  with  the  elemental  energy 
summarized  in  the  paragraphs  preceding  Eq.(3),  yielding. 


AU  = 


NplAI 

k 


yro 


['  (.P.j  +  ^^o)-F(jP  =  Pf)‘iPci 


fPf+^  1 

+  ^  {Pcj  +  ^(P/  +  1  -  M  )-F{p  =  Pf)  dp 

J  V  r-\-l—VTO  ^ 

/‘Py+l—yro 

/  (Fcj-Pf)F{p  =  Pf)dpcj 

Jo 


/p^  +  l-yro 

NplAj 
k 


‘py+1— yro 


1  rPf+^ 

/p^+l-pro 

NplAj 


F{p  =  Pf)  dpcj  +Pf  F{p  =  Pf)  dpcj 
Jo 

1  /*P/+1  1  /*P/+1 

+  9  /  ^  F{p^pf)dpcj  +  ::  ^  iPcj-PffF{p  =  Pf)dp 

^  J V t+l-yro  ^  '^P/+l“yTo 

-^[/s  —  h(p/)]  +  ^  Sj^[-^3{p/)  +  h]  +  ■^3(P/)] 


-  j  [(exp(-C^)  -  (1  -  pro)  •  ea;p(-C^(l  -  !?ro)^)]) 


(Ma) 


The  right  hand  side  of  Eq.(ll)  becomes, 
rV(p=p/) 

'V(p=0) 


fV(p=Pf)  P  rpoVTo  ^  rPf  ^  1 

/  pdV  =  PfV{p  =  Pf)  -  /  V(0  <  p  <  PoyTo)dp  +  /  ^  V{poyTo  <  Pf)dp\ 

JV(p=0)  ^*^0  J  PoVto 

PoAV 


Vto  + 1 


(yro  +  l)/-  .  T  rrr  \  ,  yTo(yTO-l)  r  ^  T 

- — t: — -iPf  +  I)  •  h{Pf)  +  - - ^ - h-—yTQ-  h 


ttA 


5l 


(146) 


/3/TO  — 1  4  2  . 

hip)  dp  +  :^ih  -  h)  -  ':^yT0  •  exp{-c  ) 

It  should  be  noted  that  the  factors,  Np^Ag/k,  in  Eq.(14a)  and  poAV/iyTo  +  1),  in 
Eq.(14b)  are  identical,  thus  dropping  out  of  the  conservation  of  energy,  Eq.(ll),  as  conunon 
factor. 

The  pkV  equation  constructed  from  the  mechanistic  model,  Eq.(13b),  contains  five 
unknowns  (A,  po,  AV,  Vu  {or  Vl),  ^ro)>  ffi©  magnitudes  of  which  are  determined  by  min¬ 
imizing  the  differences  between  Eq.(13b)  and  a  specified  data  set  based  on  the  method 
of  least  squares,  under  the  constraint  imposed  by  the  conservation  of  energy,  Eq.(ll)  and 
Eq.(14).  Because  the  p-V  equation  consists  of  two  equations,  one  for  the  high  pressure 
region  and  the  other  for  the  low  pressure  region,  and  also  because  algebraic  equations  re¬ 
sulting  from  the  application  of  the  method  of  least  squares  are  non-linear,  a  computational 
program  is  developed  to  find  the  five  unknowns.  The  program  requires  a  set  of  initial  guess 
values  for  the  five  unknowns.  The  parameters.  A,  po,  AV,  Vy,  of  the  error  function  p-V 
equation,  Eq.(lb),  are  used  for  initial  values  with  the  initial  value  for  the  fifth  unknown, 
yroi  being  set  to  zero.  The  program  employs  the  Newton-Raphson  iterative  technique 
around  the  value  of  yro  to  minimize  the  errors  between  Eq.(13b)  and  the  data  points 
while  conservation  of  energy  is  satisfied  exactly,  until  the  five  unknowns  converge  to  a  set 
of  solutions. 

Discussion  of  Results 

We  begin  with  physical  interpretations  of  parameters  of  p-V  equations  in  terms  of  the 
mechanistic  model.  The  first  equation  in  Eq.(12)  is, 

AV  =  NVoivTo  +  1) 

=  Ar(yo  +  A,yr). 

Noting  that  yro  is  a  ratio  of  the  piston  displacement  by  elastic  tissue  distension  to  that 
by  alveolar  recruitment,  AV  of  the  error  function  p-V  equation  is  the  maximum  possible 
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volume  available  for  inflation,  and  is  related  to  the  mechanistic  model  as  a  product  of  the 
total  number  of  elements,  N,  and  the  elemental  volume  available  for  inflation  through  both 
the  ‘pop-open’  mechanism,  (corresponding  to  the  volume  increase  due  to  the  alveolar 
opening),  and  the  piston  displacement,  AsVt  (corresponding  to  the  elastic  wall  distension  of 
TRS).  Vu  is  related  only  to  the  solution  of  the  high  pressure  region  as  Vi;  =  V  (p  oo). 
On  the  other  hand,  under  the  two-region  p-V  equation  of  the  mechanistic  model,  the 
definition  of  AV  needs  to  be  elaborated.  Since  Vl  7^  V  (p  -  00)  in  the  lower  pressure 
solution  of  the  mechanistic  model,  AV  should  be  interpreted  as  the  maximum  possible 
volume  change  when  the  high  pressure  solution  is  extended  into  the  low  pressure  region. 

The  second  equation,  po  =  {^/Aa)  yo,  indicates  that  the  pressme  at  the  midpoint  of  the 
p-V  cmve  is  an  equivalent  pressure  required  to  displace  the  piston  against  the  spring  force 
over  the  pop-open  displacement  of  yo.  It  may  be  rewritten  as  poVb  =  kyl;  therefore,  PqVq/2 
is  the  spring  energy  required  to  displace  the  piston  by  the  amount,  yo.  This  observation 
implies  that  the  pressure,  po,  is  related  to  both  the  alveolar  recruitment  (through  po)  and 
the  elastic  tissue  distension  (through  k).  A  higher  magnitude  of  po  implies  a  larger  value 
of  the  spring  constant  9wall  elasticity)  and/or  a  greater  amount  of  energy  required  to 
pop-open  the  elements. 

The  non-dimensional  parameter,  A,  is  related  to  the  parameter,  C,  of  the  mechanistic 
model  through  the  third  equation  in  Eq.(12),  C  —  (■v/^/4)  A.  As  may  be  seen  from  Eq.(9), 
the  parameter,  C,  appears  as  a  factor  in  the  function,  hip).  Since  the  function,  hip)i  is  a 
monotonically  increasing  function  of  p,  an  increase  in  volume,  V,  becomes  more  sensitive 
to  a  change  in  pressure  when  the  magnitude  of  A  is  larger.  The  observations  above  may 
be  further  extended  in  terms  of  the  distribution  of  elements  over  the  critical  pop-open 
pressme,  Eq.(13c).  The  number  distribution  of  elements  is  a  normal  distribution  with  its 
mean  at  p(=  p/po  -  1)  and  a  standard  deviation,  a,  which  is  proportional  to  1/A.  Since 
the  peak  of  the  distribution  is  located  at  =  p,  the  rate  of  increase  in  the  number  of 
open  elements  increases  (decreases)  for  p  <  po  (p  >  Po);  an  observation  consistent  with 
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the  fact  that  po  is  a  pressure  at  the  inflection  point  of  the  error  function  p-V  equation. 
A  larger  value  of  A  indicates  a  smaller  standard  deviation,  indicating  a  higher  peak  in 
number  density  and  a  sharper  distribution. 

The  p-V  equation,  Eq.(13b),  of  the  mechanistic  model  consists  of  the  low  pressure 
solution  in  which  all  open  elements  are  unsaturated,  and  the  high  pressure  solution  where 
some  elements  are  saturated  (fully-distended).  The  equation  has  three  pressure  -dependent 
terms.  A  term  proportional  to  p  in  the  equation  for  the  low  pressure  region  (the  third 
term)  is  due  to  the  elastic  distension  of  the  elements  that  are  open  at  p=0.  Two  other 
pressure-dependent  functions  are  originating  from  volume  changes  due  to  the  piston 

displacement,  and  which  results  from  both  the  pop-open  volume  and  the  piston 

displacement.  The  former  is  symmetric  with  respect  to  p{—  p/po  —  1)  =  0,  i.e.  h^p)  = 
hi—p)'}  while,  the  latter  is  antisymmetric,  i.e.  h{{p)  —  —h{—p)-  Furthermore,  the  p-V 
equation  in  the  high  pressure  region,  V  {po-  yxo  ^  p)i  is  independent  of  the  magnitude  of 
yro]  while,  the  solution  Vp  is  sensitive  to  the  magnitude  of  yro  in  both  the  low  and  the 
high  pressme  regions. 

Fig.  5  shows  six  representative  data  sets  of  patients  with  ARDS  as  well  as  the  cor¬ 
responding  p-V  equation,  Eq.(13b),  derived  from  the  mechanistic  model  for  the  inflation 
(I)  process.  Fig.  6  is  a  plot  of  [the  volume  predicted  by  model-based  p-V  equation  at  a 
specified  pressme]  vs  [the  corresponding  data  volume]  for  all  inflation  data  points  from 
the  twenty  one  data  sets.  Both  figmes  show  very  good  agreements  between  the  model  and 
the  clinical  data  with  R?  for  Fig.  6  being  equal  to  0.9993.  (Various  parameters  for  all 
data  sets  are  smnmarized  in  Table  1.)  The  solid  (dotted)  cmves  in  Fig.5  are  the  solution 
of  the  low  (high)  pressure  region  (i.e.  the  first  (second)  equation  in  Eq.(13b))  with  the 
composite  solution  indicating  that  the  p-V  curve  is  not  antisymmetric  with  respect  to  po. 
However,  since  the  error  minimization  is  applied  between  the  antisymmetric  error  function 
p-V  equation,  Eq.(lb),  and  the  mechanistic  model  equation,  Eq.(13b),  the  two  cmves  are 
very  close  to  each  other  in  the  low  pressme  region  of  0  <  p  <  poyro-  Ranges  of  various 
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parameters  listed  in  Table  1  are, 

A  =  1.5  -  5.5,  po  =  13  -  31  [cmH20], 

AV  =  1  -  4  [L],  VTo  =  0.289  -  0.695. 

Since  the  boundary  pressure  between  the  low-  and  the  high-  pressure  solution,  po  •  yro,  is 
low  compared  to  the  end-of  inflation  pressure  for  most  data  sets,  the  antisymmetric  high 
pressme  solution  is  applicable  over  a  major  part  of  the  data  sets  analyzed.  It  is  also  noted 
here  that  if  the  condition  of  yro  =  0  (negligible  elastic  tissue  distension)  is  imposed,  the 
solution  of  the  mechanistic  model,  consisting  solely  of  the  solution  for  the  high  pressure 
region  becomes  identical  to  the  (antisymmetric)  error  function  p-V  equation  although  the 
conservation  of  energy  is  not  satisfled  by  the  condition.  Fig.  7  is  presented  to  show  the 
magnitude  of  the  left  hand  side  of  conservation  of  energy  divided  by  poAV  as  the  abscissa, 
and  (-l)-(the  right  hand  side  of  conservation  of  energy  divided  by  poAV)  as  the  ordinate 
for  all  data  sets  when  yxo  is  set  to  zero  and  the  parameters  (A,  po,  AF,  V/,)  of  the  error 
function  p-V  equation  are  used  for  the  evaluation.  The  flgure  shows  that  the  left  and  right 
hand  side  of  conservation  of  energy  have  opposite  signs  for  all  data  sets,  indicating  that 
conservation  of  energy  is  not  satisfled  at  yTo  —  O- 

The  range  of  yro  obtained  by  the  mechanistic  model  indicates  that  the  fraction  of  total 
volume  available  for  the  pop-open  mechanism  (alveolar  recruitment),  NVq/AV,  which  is 
equal  to  1/(1  +  yro)?  ranges  between  0.59  and  0.78. 

Summary 

A  mechanistic  model  of  TRS  elements,  each  consisting  of  a  piston-spring  system,  is 
developed  to  analyze  quasi-static  pressure-volume  curves  for  the  inflation  process.  The 
model  accommodates  both  the  alveolar  recruitment  (in  terms  of  the  critical  pop-open 
pressure)  and  the  elastic  distension  of  wall  tissues  (in  terms  of  the  piston  displacement). 
Model-based  relations  (Eq.(12))  are  established  between  the  parameters  in  the  p-V  curve 
represented  by  the  error  function  equation,  Eq.(l),  and  in  the  mechanistic  model.  Under 
the  constraint  imposed  by  conservation  of  energy,  the  parameters  of  the  model-based  p-V 
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equation  is  determined  for  each  clinical  data  set  by  a  computational  minimization  proce¬ 
dure  between  the  equation  and  the  data  points,  results  of  which  show  excellent  agreements 
between  the  two  (Figs.5  and  6).  The  p-V  equation  thus  derived,  Eq.(13b),  consists  of  two 
equations;  one  for  the  low  pressure  region  where  all  open  elements  are  active  (=  unsatu¬ 
rated)  as  the  piston  of  an  element  is  yet  to  reach  its  stroke  limit,  and  the  other  for  the  high 
pressure  region  where  some  open  elements  are  saturated.  The  elemental  distribution  over 
the  critical  pop-open  pressure,  Eq.(13c),  is  a  normal  distribution  with  its  shape  (the  mean 
and  the  standard  deviation)  affected  substantially  by  the  magnitudes  of  two  parameters 
in  the  mechanistic  model,  A  and  po- 

The  present  analysis  is  for  the  inflation  process.  The  deflation  process  is  different 
from  the  preceding  inflation  because  of  the  absence  of  the  pop-open  mechanism,  and  also 
becasue  of  a  possibility  of  airway  closure.  However,  a  certain  aspect  of  the  deflation  process 
may  be  predicted  from  the  inflation  analysis;  which  will  be  discussed  in  Prat  II  as  a  validity 
test  of  the  mechanistic  model.  Wide  ranges  covered  by  the  parameters.  A,  po,  AV  and 
yro>  of  the  p-V  equation  and  the  mechanistic  model  need  to  be  interpreted  in  terms  of  the 
shape  and  the  range  of  the  p-V  curves  as  well  as  in  terms  of  the  elemental  distribution  and 
its  changes  along  the  corresponding  p-V  curve;  which  will  also  be  discussed  in  Part  II. 
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Figure/Table  Captions 

Table  1.  Summary  of  Inflation  Data  Analyses. 

Fig.l.  A  typical  quasi-static  pulmonary  pressure- volume  curve. 

Fig.2.  Error-function  p-V  equation  and  inflation  data  points. 

Fig.3.  Error-function  p-V  equation  and  deflation  data  points. 

Fig.4.  A  schematic  diagram  of  mechanistic  model  of  TRS  element. 

Fig.5.  Model-based  p-V  equation,  Eq. (13b),  vs  data  points  for  inflation  process, 
solid:  solution  for  low  pressure  region, 
dotted;  solution  for  high  pressure  region. 

Fig.6.  V  (voliune  predicted  by  model-based  p-V  equation)  vs  V  (volume  of  data) 
for  a  specifled  pressme. 

Fig.7.  AU/po  •  AV  vs  (-1)  •  /  pdVfpo  •  AV  when  yro  =  0 
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Table  1.  Summary  of  Inflation  Data  Analysis. 


Data  No 

A 

Po 

[cmH20] 

AR 

(L] 

Vl 

[L] 

Vto 

A.l. 

2.9652 

22.398 

2.3726 

-0.0709 

2. 

2.9578 

22.411 

2.3772 

-0.0723 

0.347 

B.l. 

2.7173 

21.981 

1.5612 

-0.0762 

2. 

2.7304 

21.999 

1.5559 

-0.0728 

0.359 

C.l. 

3.3532 

25.082 

1.6193 

-0.0365 

2. 

3.3664 

25.073 

1.6139 

-0.0342 

0.329 

D.l. 

1.6273 

13.324 

3.1567 

-0.5727 

2. 

1.9257 

13.999 

2.8392 

-0.3542 

0.379 

E.l. 

2.7160 

30.361 

1.6216 

-0.0768 

2. 

2.6497 

30.817 

1.6847 

-0.0923 

0.480 

F.l. 

2.6029 

23.880 

1.5066 

-0.0962 

2. 

2.6423 

23.863 

1.4874 

-0.0873 

0.365 

G.l. 

1.7288 

20.156 

3.0989 

-0.3645 

2. 

1.9277 

20.731 

2.8704 

-0.2138 

0.447 

H.l. 

1.8901 

14.959 

1.7905 

-0.2185 

2. 

2.0421 

15.405 

1.7129 

-0.1563 

0.431 

I.  1. 

3.5379 

25.248 

2.7570 

-0.0350 

2. 

3.5449 

25.232 

2.7508 

-0.0333 

0.368 

J.  1. 

2.7981 

26.208 

3.7326 

-0.1583 

2. 

2.7364 

26.830 

3.9129 

-0.1887 

.  0.474 

K.l. 

2.4296 

17.895 

1.3424 

-0.0997 

2. 

2.4708 

17.951 

1.3304 

-0.0916 

0.440 

L.l. 

1.2500 

11.592 

1.2470 

-0.2829 

2. 

1.5318 

13.213 

1.1256 

-0.1678 

0.695 

MI. 

3.2725 

29.865 

3.9075 

-0.0861 

2. 

2.9972 

30.327 

4.2463 

-0.2179 

0.626 

N.l. 

2.7915 

15.310 

1.6284 

-0.0705 

2. 

2.8046 

15.297 

1.6219 

-0.0681 

0.358 

0.1. 

1.9487 

18.374 

1.8797 

-0.2423 

2. 

2.1412 

18.837 

1.7753 

-0.1696 

0.389 

Rl. 

2.4566 

26.982 

1.3277 

-0.0787 

2. 

2.4553 

27.032 

1.3278 

-0.0773 

0.467 

Q.l. 

2.1027 

19.314 

1.3316 

-0.1279 

2. 

2.2041 

19.583 

1.2941 

-0.1004 

0.389 

R.I. 

1.1672 

13.986 

2.0685 

-0.5149 

2. 

1.6209 

16.352 

1.7396 

-0.2530 

0.406 

S.l. 

3.1381 

26.802 

3.1306 

-0.0851 

2. 

3.0894 

26.925 

3.1878 

-0.1030 

,  0.407 

T.l. 

5.4709 

30.038 

1.7695 

-0.0097 

2. 

5.4708 

30.037 

1.7694 

-0.0097 

0.289 

U.l. 

3.2818 

24.439 

2.8956 

-0.2075 

2. 

3.2819 

24.349 

2.8956 

-0.2075 

0.327 

1.  Obtained  by  applying  the  method  of  least  squares 
along  with  error  function  p-V  equation. 

2.  Results  from  the  mechanistic  model. 
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Fig. 2.  Error-function  p-V  equation  and  inflation  data  points. 
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Flg.3.  Error-function  p-V  equation  and  inflation  data  points. 


Fig.5.  Model  based  p-V  equation,  Eq.  (13b),  vs  data  points  for  inflation  process. 
Solid:  solution  for  low-pressure  region, 

Dotted:  solution  for  high-pressure  region. 
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Fig.6.  V  (volume  predicted  by  model-based  p-V  equation)  vs.  V(Volume  of  data)  for  a  specified  pressure. 
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Chapter  3.  A  Mechanistic  Model  for  Quasi-Static  Pulmonary  Pressure- Volume 
Curves:  Examination  of  Clinical  Data 

Abstract 

A  p-V  equation  is  developed  in  Part  I  based  on  a  mechanistic  model  of  a  total  respiratory 
system.  In  Part  II,  twenty  one  p-V  data  sets  of  patients  with  acute  respiratory  distress 
syndrome  are  examined  using  the  mechanistic  model,  relating  the  quasi-static  putoonary 
p-V  curve  to  the  corresponding  respiratory  conditions  in  terms  of  a  volume  increase  due 
to  alveolar  recruitment  and  due  to  elastic  tissue  distension,  the  elemental  distribution 
ranging  from  the  closed  elements  to  the  saturated  (open  and  fully-distended)  elements 
and  its  changes  with  pressure.  The  compliance  (local  gradient)  of  p-V  curves  is  shown 
to  represent  the  change  in  the  total  volume  of  saturated  elements;  while  the  pressure  at 
the  maximum  compliance  is  identified  as  the  location  where  a  maximum  rate  of  increase 
occurs  both  in  the  volume  increase  due  to  alveolar  recruitment  and  in  the  volume  increase 
due  to  an  increase  in  the  saturated  elements.  Validity  of  the  model  is  provided  by  its 
predictions  of  the  upper  volume  asymptote  and  the  maximum  possible  volume  change  of 
the  corresponding  deflation  process  which  agree  well  with  the  clinical  data. 
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Nomenclature 


>1. 

B 

C 

/,  F 

h  (i  =  1  -  5) 

k 

N 

^open 

Ngat 

^unsat 

Nj 

P 

P 

Pcj 

Pcj 

Pf 

Po 

PID 

V 


Vpop—open 


'sat 


'^open—sat  (P  — 


v„ 


V, 


L(U) 


piston  surface  area  on  which  pressure  is  acting 
{k/As)yT  =  Po  •  Vto 
—  ■y/7rA/4 

distribution  functions  (Eq.(3)) 

functions  defined  in  Eq.(2) 

spring  constant  [N/m] 

total  number  of  TRS  elements 

total  nmnber  of  open  elements 

total  number  of  saturated  elements 

total  number  of  unsaturated  elements 

number  of  elements  at  energy  level  j 

pressure  (interplemal  pressme  difference) 

non-dimensional  pressure,  p/po  —  1 

critical  pressure  at  which  an  element,  j,  ‘pops  open’. 

Pcj/PO 

pressure  at  the  end  of  inflation 

pressure  at  the  inflection  point  in  model  equation 

pressure  at  the  intersect  of  inflation  and  deflation  processes 

volume 

total  ‘pop-open  volrnne’  ’ 

total  volume  of  satmated  elements 

0)  total  volume  of  elements  open  at  p  =  0  when  they  are  all  saturated, 
volume  change  from  the  state  of  p  =  0 
lower  (upper)  boimd  of  volume 

an  (imaginary)  upper  boimd  of  volume  for  the  deflation  process 
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V 


non-dimensional  volume,  {V  —  {{Vu  +  Vl)/2)/(AV/2),  (Eq.(lb)) 


volume  of  an  element  j 

AV 

Vu  —Vl  =  NVoivTo  +  1) 

K, 

‘pop-open’  voliune  (=  A^yo) 

Vj 

piston  displacement  of  an  element  j 

yo 

It 

Co 

Vt 

piston  stroke  limit 

Vto 

Vt/vo 

Greek  symbols: 

Q 

constant  of  proportionality 

A 

apo A V(non-dimensional  parameter)  (Eq.(lb)) 

a 

(8/7r)^/^/A,  Standard  deviation  (Eq.(3b)) 

Od 

(7  •  po 

Superscript; 

related  to  a  single  TRS  element 

d 

deflation  process 

Acronym: 

ARDS 

acute  respiratory  distress  syndrome 

TRS 

total  respiratory  system 
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Introduction 

In  Part  I,  the  error  function  p-V  model  equation  is  shown  to  agree  well  with  clinical 
p-V  data.  The  equation  and  the  corresponding  non-dimensional  form  are, 

where  AV  —  Vu  —  VliVu  =  upper  volume  asymptote,  Vl,  =  lower  volume  asymptote, 

Po  =  a,  pressmre  at  the  inflection  point  of  the  curve,  C  =  y/nKfA, 

V  =  [V-{Vu+  Vl)/2]/{AV/2),  p  =  (p/po)  -  1. 

The  mechanistic  model  of  a  TRS  element  developed  in  Part  I  is  a  piston-spring- 
cyhnder  system  with  the  alveolar  recruitment  and  the  elastic  tissue  distension  represented 
respectively  by  the  critical  pop-open  pressure,  Pcj,  and  by  the  displacement  of  piston 
against  the  spring  force,  pj  =  Ag(p-pcj)/k  (Ag  =  piston  surface  area,  k  =  spring  constant). 

Based  on  the  error  function  p-V  equation  and  the  mechanistic  model,  and  allowing 
for  a  distribution  of  elements  over  pcj,  the  following  model-based  p-V  equations  as  well  as 
the  corresponding  distribution  function  of  TRS  elements  are  derived. 


V  (  0  <  p  <  Po  •  Pto  )  = 


Vu+Vl 


Pro  + 1 


im  +  1  -  M 


,  (2a) 


T..  ^  ^  X  Vu  +  Vl  .  AV  , 

V  (po  •  Pro  <  p)  = - 2 - ~ 


where 

Ji=er/(C'),  l2{p)  =  exp{-C^p^)-exp{-C^),  l3{p)  =  erf{Cp), 

h  =  erf{C(l-yTo)),  h  =  exp{-C^  {yro  -  if)  -  exp{- C"^), 

CD  =  (8/7r)5  po/A,  a  =  (S/tt)  Pcj  =  PcjlPo- 
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Following  the  development  of  the  mechanistic  TRS  model  and  the  discussion  on  relations 
between  the  error  function  p-V  equation  and  the  the  model-based  p-V  equation  as  well  as 
between  parameters  of  the  two  equations  in  Part  I,  our  discussions  here  concentrate  on 
various  results  that  may  be  derived  from  applications  of  the  model  to  the  clinical  p-V  data 
sets. 

Relationship  between  Inflation  and  Deflation  Processes 

Although  this  report  is  focused  on  the  mechanistic  model  for  the  inflation  process, 
there  exist  certain  relations  between  the  inflation  and  the  deflation  process  that  may  be 
evaluated  from  the  present  inflation  analyses.  We  consider  a  general  case  in  which  a 
quasi-static  inflation  process  proceeds  to  a  pressure,  pio  (=  end-of-inflation  pressure  = 
initial  pressure  of  the  corresponding  deflation  process) ,  followed  by  a  quasi-static  deflation 
process.  In  terms  of  the  mechanistic  model,  TRS  elements  at  p  =  pm  with  its  critical 
pop-open  pressure  less  than  zero  {pcj  <  0)  are  still  closed  and  have  not  contributed  to 
the  volume  change  dining  the  inflation  process  from  p  =  0  to  p  =  p/£j;  hence,  we  may 
postulate  that  only  those  elements  that  are  open  at  p  =  p/£)  participate  in  the  deflation 
process  to  follow.  Therefore,  Vy  (=an  (imaginary)  upper  bound  of  volume  for  the  deflation 
I>-V  curve)  may  be  viewed  as  the  volume  which  would  be  attained  if  the  elements  that  are 
open  at  the  end  of  the  inflation  process,  p  =  P/d>  were  all  fully  saturated;  i.e. 

Wo(l  +  Vto)  •  =  Pid)  dPcj  -  Vp{p  =  Pjd) 

The  first  term  on  the  right  hand  side  is  the  inflation  volume  at  p  =  p/o.  The  second  integral 
term  is  the  volume  summed  over  all  open  elements  at  p  =  piD  when  they  are  saturated, 
and  the  last  term  is  the  actual  volume  increase  in  the  inflation  process  from  p  =  0  to 
P  =  Pi D  with  the  two  terms  in  the  square  bracket  together  representing  a  volume  increase 
above  V(p  =  p/p)  if  all  open  elements  at  p/p  were  saturated.  Under  the  assumption 
that  the  magnitude  of  p/p  is  greater  than  B  (=  PoVto)  which  is  valid  for  all  data  sets 
analyzed,  Eqs.(2,3)  along  with  Eq.(9)  of  Part  I  for  Vp  are  used  to  evaluate  the  right  hand 
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VS  =  V(p  =  p,d)  + 


f 


side,  yielding 


7rA(H-yroj  2 

AV 

+  [1  -  (2  -  yro)  • /i  +  (1  -  ^To)  •  i4]  (4) 

The  clinical  data  sets  made  available  to  us  contain  both  the  inflation  and  the  deflation  p-V 
cmves  for  each  patient  with  ARDS;  however,  the  p-V  curves  are  obtained  separately  for  the 
inflation  and  the  deflation  process.  (See  [1]  for  the  procedure  of  data  acquisition.)  Fig.l 
shows  inflation  (unfilled)  and  deflation  (filled)  data  points,  as  well  as  the  corresponding 
inflation  (I)  and  deflation  (D)  curves  for  a  typical  data  set  we  examined.  The  inflation 
curve  in  Fig.l  is  Eq.(2)  of  the  mechanistic  model;  while,  the  deflation  curve  is  obtained 
by  straight  applications  of  the  method  of  least  squares  between  data  points  and  the  error 
function  p-V  equation,  Eq.(l).  As  may  be  observed  from  Fig.l,  the  end-of-inflation  point  is 
quite  difierent  from  the  initial  deflation  point  for  most  data  sets.  To  acconunodate  the  data 
into  the  analysis  based  on  Eq.(4),  the  initial  deflation  data  point  is  translated  horizontally 
imtil  it  meets  the  inflation  p-V  curve,  the  pressure  value  of  which  is  then  defined  as  pjd 
in  Eq.(4),  as  indicated  in  Fig.l,  implying  that  the  deflation  curve  preceded  by  an  inflation 
curve  is  assumed  to  be  the  same  as  the  deflation  curve  of  data  sets  horizontally  translated 
until  the  beginning-of-deflation  data  point  is  on  the  inflation  curve. 

Fig.  2  presents  of  Eq.(4),  predicted  from  the  mechanistic  model  of  the  inflation 

process,  plotted  against  of  the  error-function  p-V  equation,  Eq.(l).  (For  a  complete 
list  of  numerical  results  relevant  to  the  analysis,  see  Table  1.)  A  maximum  and  a  minimum 
of  a  difiference,  (jE'g.(4))  —  {Eq.{l)),  are  0.1113  [L]  and  -0.0352  [L]  respectively  with 

an  average  of  the  difference  =  0.0460  [L].  Agreements  are  very  good  in  view  of  the  fact 
that  Eq.(4)  predicts  the  upper  volume  asymptote  of  the  deflation  process  in  terms  of  the 
conditions  predicted  by  the  mechanistic  model  of  the  corresponding  inflation  process;  thus 
indirectly  supporting  a  certain  degree  of  validity  of  the  mechanistic  model.  Also,  the  fact 
that  the  magnitude  of  pjD  is  determined  from  the  horizontal  translation  of  the  deflation 
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curve  indicates  that  the  deflation  process  may  be  relatively  insensitive  to  the  inflation 
history  prior  to  p  =  p/o.  On  the  other  hand,  Eq.(4)  for  contains  piD  Sts  variable  for 
specified  inflation  conditions,  indicating  that  the  magnitude  of  (i.e.  the  shape  of  the 
deflation  curve)  changes  as  the  end-of-inflation  pressure  is  varied.  A  similar  statement 
has  been  made  previously  by  Jonson  [2].  A  closer  examination  of  Fig.  2  indicates  that 
the  mechanistic  model  slightly  underpredicts  compared  to  that  of  the  error-function 
equation  for  most  of  data  sets.  This  could  indicate  either  a  quantitative  limitation  of  the 
mechanistic  model  or  the  effects  of  the  inflation  process  preceding  the  deflation.  If  yxo  is 
set  to  zero,  Eq.(4)  is  reduced  to 


VS  iVTo  =  0)  =  Vu  -  AV  (/i  -  /3(p/o))/2.  (4a) 


Fig.3  plots  Vy  of  Eq.(4a)  vs  of  the  error  fimction  p-V  equation  for  the  deflation  process. 
Agreements  are  fairly  good  between  the  two.  Results  presented  in  Fig.3  reflects  that  the 
p-V  curve  is  relatively  insensitive  to  the  magnitude  of  yro]  a  reason  why  the  antisymmetric 
error  function  p-V  equation  (for  which  pro  is  zero)  fits  well  with  p-V  curves. 

The  error  function  p-V  equation  fits  well  not  only  with  the  inflation  but  also  with  the 
corresponding  deflation  processes  as  shown  in  Part  I.  Therefore,  regardless  of  the  actual  de¬ 
flation  process,  we  may  define  AV^  (A V  of  the  deflation  process)  as  the  maximum  possible 
volume  change  of  a  specified  TRS  during  the  deflation  process;  which,  in  our  mechanistic 
model,  yields  the  following  equation  for  AV*^  in  terms  of  the  inflation  parameters: 


poo 

AV<^  =  VoiyTo  +  l)  /  N-F{p  =  pjD)dpcj 
Jo 


AV 


(1  -I-  hiPin))- 


(5) 


Fig.4  is  AV  of  deflation,  AV*^,  predicted  by  Eq.(5)  plotted  against  the  corresponding  AV^ 
of  error-function  p-V  equation  determined  by  the  method  of  least  squares.  (See  Table  1 
for  numerical  values.)  Agreements  between  the  two  are  reasonably  good  for  a  majority 
of  data  sets,  except  for  six  data  sets  shown  in  filled  circles  accompanied  by  alphabetical 
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data  numbers.  Shown  in  Fig.5  for  the  six  data  sets  are  their  deflation  data  points  and 
two  p-V  curves;  one  (dotted)  for  the  error  function  p-V  equation  with  the  method  of  least 
squares  applied  to  determine  (A"^,  Po>  V^),  and  the  other  (solid)  for  the  error  function 
p-V  equation  with  (AV*^,  V^)  determined  from  Eqs. (4,5)  and  (A*^,  pg)  determined  by  the 
method  of  least  squares.  Two  curves  are  diflferent  in  their  approaches  to  different  high  and 
low  asymptotes.  Although  agreements  of  the  solid  curves  (with  two  adjusting  parameters) 
with  the  data  points  are  not  as  good  as  that  of  the  dotted  curves  (with  four  adjusting 
parameters),  the  errors  are  small  for  the  solid  curves  in  view  of  the  fact  that  the  magnitude 
of  AV*^  is  quite  different  between  the  two  curves,  indicating  that  a  better  understanding  of 
relations  between  the  p-V  equation  and  the  corresponding  intra-respiratory  changes  helps 
interpret  various  characterisitcs  of  p-V  curves  accurately. 

Interpretation  of  Inflation  i>-V  Curves  based  on  Mechanistic  Model 

Fig.  6  shows  ranges  covered  by  all  data  sets  analyzed  in  terms  of  po  (the  inflection 
pressure  of  the  high-pressure  solution),  AV  (maximum  volume  available  for  inflation  in 
the  high-pressure  solution)  and  l/(l-l-pro)  (the  fraction  of  total  elemental  volume  available 
for  the  pop-open  mechanism  (alveolar  recruitment),  NVq/AV),  all  plotted  against  the  non- 
dimensional  parameter,  A.  The  data  sets  with  their  alphabetical  data  numbers  indicated 
in  the  figure  are  those  to  be  analyzed  in  detail  in  comparative  analyses  to  follow.  (Various 
parameters  of  the  six  data  sets  are  reproduced  as  Table  2.  Parameters  of  all  data  sets 
are  listed  in  Table  1  of  Part  I.)  The  range  of  A  is  between  1.5  and  3.6  except  for  Data 
T  (A  =  5.47).  The  six  data  sets  (B,  E,  M,  N,  R,  T)  are  different  from  each  other  in  the 
following  ways: 

1.  Data  set  B,  E  and  N  have  roughly  the  same  magnitude  in  both  A  (=~  2.65  —  2.80)  and 
AV  (=~  1.55  —  1.68)  with  significantly  different  values  for  po. 

2.  Data  set  E  and  T  have  substantially  different  values  of  A  with  po  and  AV  being 
approximately  the  same  in  magnitude. 
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3.  N  and  R  also  show  similar  characteristics  of  being  different  in  A  eind  common  in  pq  and 
AV";  however,  the  pair  covers  lower  range  in  both  A  and  po  compared  to  E  and  T. 

4.  E  and  M  are  different  in  terms  of  the  magnitude  of  AV  with  pqe  ~  Po  Mj  ~  A  a/. 

5.  T  and  R  represent  the  data  sets  with  very  high  and  low  values  of  A,  respectively. 

Figs.  7  and  8  are  p-V  and  the  corresponding  non-dimensional  p  —  V  curves  over  their  ranges 
of  measurement  for  Data  Set  B,  E  and  N  in  (a),  Data  Set  E  and  T  in  (b),  Data  Set  N 
and  R  in  (c)  and  Data  Set  E  and  M  in  (d). 

Referring  to  Fig.  7  and  noting  that  the  dotted  (broken)  vertical  lines  are  the  location 
of  p  =  Po  (po  •  yTo))  it  be  observed  that  the  range  for  the  low  pressure  solution, 
0  <  p  <  Po  •  yrO)  in  which  all  elements  are  active  (unsaturated),  is  very  narrow  compared 
to  the  range  for  the  high  pressure  solution  with  an  exception  of  Data  set  M,  for  which  the 
measurement  does  not  reach  the  inflection  pressure,  po,  with  po  •  yro  >  15  [cmH20]. 

The  corresponding  non-dimensional  (p  —  V)  curves,  based  on  the  definitions  of  V  and  p  in 
Eq.(l),  represent  Eqs.(2a,b)  in  the  following  normalized  form: 


( -1  <  P  <  yro  -  1 )  = 

,  2  ( T  i-\  ,  (l--^i)/-  ,1  -  \  yro  - 1  r  ,  2(75  -72(p))\ 

+  +  /4+— - )’ 

V  (yro  —  1  <  y) 


=  h{p), 


(6) 


The  normalization  of  volume  transforms  the  two  volume  asymptotes,  Vu  and  Vl  into 
V  =  +1  and  V  =  —  1  respectively;  while,  the  pressure,  p  =  po,  at  the  inflection  point 
is  transformed  into  p  =  0.  With  both  the  location  of  po  and  the  volume  asymptotes 
made  common  to  all  p-V  curves,  the  resulting  non-dimensional  representations  in  Fig.8 
are  characterized  by  a  single  non-dimensional  parameter,  A.  The  parameter.  A,  is  twice 
the  maximum  local  compliance  at  p  =  po  {dV/dp{p  =  0)  =  A/2).  Since  the  compliance 
is  maximum  at  the  origin  of  p  —  V  diagram,  the  first  quadrant  (V,  p  >  0)  in  Fig.8  is  a 
region  of  decreasing  local  compliance  with  pressure;  while,  the  third  quadrant  (V,  p  <  0) 
is  a  region  of  increasing  local  compliance  with  pressure.  The  origin  (p  =  0,  V  =  0)  of 
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dimensional  p-V  curves  is  transformed  into  {p  —  —1,  V{V  =  0))  on  a  p-V  curve;  hence, 
the  physiological  lower  limit  of  p  is  —1. 

Fig.8  (a)  compares  the  three  data  sets,  B,  E  and  N,  among  which  the  magnitude  of  po  is 
substantially  different  with  A  and  AV  being  approximately  the  same  in  magnitude.  Three 
curves  are  very  close  to  each  other  because  values  of  A  are  similar,  and  the  difference 
between  the  three  appears  as  the  extent  to  which  the  p-V  cmves  are  measured  in  the 
region  of  decreasing  compliance  with  pressure.  In  Fig.8  (b)  and  (c)  differences  between 
the  two  data  sets  occur  in  the  magnitude  of  A  (At  >  A^,  A;v  >  A/j),  resulting  in  the  T- 
and  N-  cmves  above  the  E-  and  R-  cmves  respectively  in  the  first  quadrant.  Since  the 
magnitude  of  po  for  the  data  set  T  and  E  are  very  high  compared  to  those  for  N  and  R, 
the  region  of  decreasing  compliance  covered  by  the  T-  and  El-  cmves  are  narrower  than  N 
and  R.  Because  both  po  and  AV  are  similar  in  magnitudes  between  the  two  data  sets  in 
Fig.8(b)  and  (c)  the  shape  ofp— .V  cmves  is  very  similar  to  the  corresponding  p-V  cmves. 
In  Fig.8  (d)  two  data  sets  with  a  high  value  of  po  (~  30  [cmH20])  are  shown.  For  the  data 
set  M  the  high  value  of  po  combined  with  a  high  value  of  AV  limit  the  measmed  range 
of  the  p-V  cmve  in  the  region  of  increasing  compliance  only,  compared  to  the  data  set  E 
with  a  smaller  magnitude  for  AV. 

Although  the  p  —  V  diagram  helps  distinguish  differences  among  p-V  cmves  and 
effectively  bring  out  various  characteristics  of  each  p-V  cmve,  it  is  the  information  from 
the  elemental  distribution  that  relates  various  parameters  of  p-V  cmves  to  TRS  conditions 
quantitatively.  On  the  normalized  p  —  V  plane  of  Fig.8  the  local  compliance  at  p  =  0 
(dV /dp  (p  =  0))  increases  with  A;  while,  as  the  standard  deviation,  cr,  is  proportional  to  1/ A 
in  the  normalized  number  distribution,  Eq.(3b),  the  distribution  becomes  sharper  and  has 
a  higher  peak  as  A  is  increased.  Fig.9  is  a  plot  of  the  number  distribution  (not  normalized) 
vs  the  critical  pop-open  pressme,  Pcj  in  [cmH20],  for  the  six  data  sets  analyzed  in  Figs. 7 
and  8.  The  number  distribution,  dNj /N  •  dpcj  >  is  a  fraction  in  the  number  of  elements,  the 
critical  pop-open  pressme  of  which  ranges  between  pcj  and  Pc  j  +  dpc  j  in  [l/cmH20].  The 
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corresponding  equation  is  Eq.(3a),  which  indicates  that  the  maximum  number  of  elements 
are  present  at  Pcj  =  p  —  Po  with  its  magnitude  equal  to  1/v^  •  cd  (=  A/4  •  po).  As  ctd 
decreases  the  distribution  becomes  sharper  (oc  1/ctd)  and  its  peak  value  (a  ao)  larger. 
(See  Table  2  for  the  magnitude  of  crjr>  for  each  data  set.)  Two  distributions  are  shown  in 
Fig.9  for  each  data  set,  one  at  p  =  po  (with  its  peak  at  p  =  0)  and  the  other  at  p  =  p/  (end- 
of-inflation  pressure)  (with  its  peak  at  p  =  p/  —  po)-  The  vertical  line,  p  =  po  (p/)  for  the 
distribution  at  p  =  po  (p/)  indicates  the  pressure  above  which  the  distribution  is  truncated. 
The  dotted  parts  of  the  curves,  as  discussed  in  Part  I,  correspond  to  the  elements  that 
are  open  at  p  =  0.  The  normal  distribution  truncated  at  Pcj  =0  and  Pcj  —  p  (i.e.  the 
solid  part  of  the  curves  in  Fig.9)  shifts  to  the  right  with  an  increase  in  p  as  more  elements 
become  open.  When  p  <  po,  the  peak  of  the  distribution  lies  in  the  negative  range  of  pcj. 
It  should  be  noted  that  in  Data  set  M  (Fig.9  (c))  the  distribution  at  p  =  p/  lies  below  that 
at  p  =  Po  because  the  measured  range  never  reached  p  =  po- 

An  integral  of  the  distribution  function  over  the  critical  pop-open  pressure  in  Fig.9  should 
yield  various  fractions  in  number  of  elements  (depending  on  the  upper  and  lower  limits  of 
the  integral)  at  each  quasi-static  state.  Also  an  integral  of  a  product  of  the  distribution 
fimction  and  the  elemental  volume  over  the  critical  pop-open  pressure  should  provide  us 
with  such  quantities  as  the  volume  change  due  to  alveolar  recruitment,  due  to  elastic  wall 
distension  and  due  to  an  increase  in  the  saturated  elements.  The  following  equations  may 
be  obtained  for  changes  in  the  number  fractions: 

Fraction  of  the  number  of  open  elements  at  p  =  p:  Nopenip  —  p)/N 

F(p)dp,i 

Fraction  of  the  number  of  saturated  elements  at  p  =  p:  Ngat.iP  —  P)/^ 

N.at.(p  =  P)  (_  \  /  [/4  +  /3(p)1/2,  if  W  •  W)  <  P. 

N  y‘Jo  PWV^>j-[o,  ifO<p<po-!?ro. 


j  =[/i+/3®)/2,  far  p>0,  (7a) 
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Fraction  of  the  number  of  unsaturated  elements  at  p  =  p:  Nunsat.{P  =  P)/^ 

^unsat.ip  ~  P')  f _  ^open  ^sat.  ^  _  f  \Il  if  Po  '  VTO  ^  P> 

N  y  N  ^  \  [/i  +  /3(p)]/2,  if  0  <  p  <  po  •  VTO- 

Fraction  of  the  number  of  open  elements  at  p  =  0;  Nopen{p  =  0)/N 

Nopenjp  =  0) 

It  should  be  noted  that  (1.)  the  equation  for  Nopen{p)/N  does  not  include  the  elements 
that  are  open  at  p  =  0.  (2.)  the  equation  for  Naat.ip)/^  does  not  account  for  the  elements 
that  are  open  at  p  =  0  and  saturated  subsequently,  (3.)  Nopen{p  =  0)/-^  is  a  function  of 
a  single  parameter,  A.  (4.)  the  number  fraction  of  unsaturated  elements,  Nunsat.ip)/N,  is 
independent  of  pressure  for  po  •  Pro  <  P,  indicating  that  as  more  elements  are  recruited  in 
the  region  of  the  high  pressure  solution  the  same  number  of  elements  are  saturated. 
Similarly,  the  following  equations  are  for  volume  changes  as  pressure  is  varied: 

Total  pop-open  volume:  Vpop~open{p  =  p) 

v„-„pe„(p  =  p)  NVo  ■  ^  +  ^3®]  ^  o-  (S") 


F{p=-l)dpcj  =|l-/il/2 


Total  volume  of  saturated  elements:  Vsat.ip  =  p) 

V.a..(P  =  P)  ^  s  W„(l  +  Jro)  ■  j 

^  f  Al/  •  [/4  -I-  /3(p)]/2,  if  Popro  <  P,  /g^'v 

\  0,  if  0  <  p  <  Po  •  Pto-  ^ 

Total  volume  of  elements  open  at  p  =  0  when  they  are  all  saturated:  Vopen-sat.ip  =  0) 

K,.„-.^.(p  =  0)  f  s  Worn  •  ^ 

^  |T£-(1  -  Ji)  for  poyTo<P-  (8c) 

2  1  +  Pto 

A  number  fraction  of  open  elements.  Nopen  total  (p  =  p)/N,  may  be  derived  from  Eqs.(7a) 


and  (7d)  as. 

Nopen  total  (P  ~  P)  f _  NopeniP  ~  0)  Nppen (jP  P)  A  _  ^  -^sCp) 

N  V"  ^  ^  /  2 


(9) 
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The  function,  hip)  defined  in  Eq.(2),  is  an  error  function,  the  magnitude  of  which  depends 
on  A  and  po-  Therefore,  the  number  fraction  of  all  open  elements.  Nopen  total  (p  =  p)/N, 
has  the  point  of  antisymmetry  at  p  =  po  with  two  asymptotes  of  0  and  1.  Fig.  10  shows 
Nopen  total  (p  =  p)/N  VS  p  of  the  six  representative  data  sets  in  solid  lines.  Due  to  the 
antisymmetry  with  respect  to  p  =  po,  the  firaction.  Nopen  total  (p  —  p)/N,  is  exactly  equal 
to  0.5  when  p  is  equal  to  po  as  the  peak  of  the  distribution  is  located  at  p  =  0.  (See 
Fig.9.)  Since  the  location  of  p  =  poPro  is  the  boundary  between  the  low-  and  high- 
pressure  solutions  of  the  p-V  equation,  all  open  elements  are  still  active  and  unsaturated 
for  p  <  PoPto-  The  fi'action  of  open  elements  at  p  =  PoPTO  is  less  than  0.2  for  four 
data  sets  other  than  Data  sets  R  and  M  for  which  the  firaction  is  ~  0.25.  As  pressme 
increases  beyond  poPTO)  some  of  the  open  elements  begin  to  be  saturated,  the  fraction 
of  which,  NsatXP  =  P)/N  of  Eq.(7b)  is  plotted  in  Fig.  10  in  broken  lines.  The  number 
firaction  saturated  depends  on  pro  in  addition  to  A  and  po-  The  difference  between  the 
two  fractions  plotted  in  Fig.  10  varies  with  the  magnitudes  of  the  three  parameters  with 
{Nopen  total /N  —  NgaiJN)  at  a  Specified  pressure  ranging  from  ~  0.28  for  Data  R  to  less 
than  0.01  for  Data  T.  It  should  also  be  noted  that  the  two  curves  are  parallel,  indicating 
that  the  rate  of  increase  in  the  nmnber  of  saturated  elements  is  equal  to  the  corresponding 
rate  of  opening  elements,  once  presstue  exceeds  PoPto- 

The  niunber  firaction  of  open  elements  at  p  =  0,  Nopen{p  =  0)/-^5  of  Eq.(7d)  rep¬ 
resents  the  percentage  of  elements  that  only  experience  elastic  displacement,  which  may 
be  interpreted  as  the  elemental  fraction  representing  a  non-alveolar  part  of  TRS  such  as 
airway  tissues  as  well  as  a  dysfunctional  alveolar  part  which  does  not  respond  to  the  re¬ 
cruitment.  Fig.  11  depicts  the  firaction  as  a  function  of  A,  the  only  parameter  affecting  the 
firaction.  As  Eq.(7d)  shows,  the  firaction,  NopenijP  =  0)/-^)  has  two  asymptotes  of  0  (as 
A  — oo)  and  0.5  (as  A  — >  0).  Fig.ll  indicates  that  the  number  firaction  of  open  elements 
at  p=0  is  very  sensitive  to  the  magnitude  of  A  as  its  value  drops  firom  0.5  to  ~  0.1  when  A 
is  changed  firom  0  to  2.  It  should  also  be  mentioned  here  that  the  total  volume  of  elements 
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open  at  p  =  0  when  they  are  all  satiu-ated,  Vopen-sat.{P  =  0)  of  Eq.(8c),  depends  on  three 
parameters,  A,  AV  and  yro,  of  the  mechanistic  model.  The  last  column  of  Table  2  lists 
Vopen-sat.ip  =  0)  for  the  six  data  sets.  Although  their  magnitudes  are  small  (0.08  [L]  or 
less),  the  value  varies  substantially  among  the  data  sets. 

As  mentioned  in  Part  I,  the  optimization  of  the  ventilator  strategy  is  required  for 
patients  with  acute  lung  injury  in  intensive  care  units  in  terms  of  pressure-  and  volume- 
ranges  [3  —  5],  taking  into  account  such  considerations  as  a  change  in  the  recruited  volume 
with  pressmre,  a  rate  of  recruitment  and  overdistension  of  the  respiratory  tissues.  Results 
of  the  mechanistic  model  analyses  relevant  to  the  respiratory  ventilation  are  presented  in 
Figs.  12  and  13.  The  volume,  Vpop-open  of  Eq.(8a),  represents  a  volume  change  of  TRS 
due  to  the  pop-open  mechanism  (alveolar  recruitment)  only  (that  is,  excluding  the  volume 
change  due  to  the  displacement  of  piston  (elastic  tissue  distention));  on  the  other  hand, 
the  volmne,  Vsat.  of  Eq.(8b),  is  the  total  volume  of  satmated  elements  of  TRS.  They 
both  increase  as  pressure  increases  along  an  inflation  path.  Sketched  in  Fig.  12  are  p  —  V 
curves  (solid,  Eq.(2)),  p  -  Vpop-open  curves  (dotted,  Eq.(8a))  and  p  -  Vsat  curves  (broken, 
Eq.(8b))  for  the  six  data  sets.  The  vertical  broken  line  in  the  figure  is  the  location  of 
p  —  po.  The  intersect  between  the  p  —  Vsat  curve  and  the  x  axis  in  the  figure  is  the  pressure 
at  the  boundary  (=  PoVto)  between  the  low  pressure  p-V  solution  (in  which  all  open 
elements  are  active  (unsaturated))  and  the  high  pressure  solution  (in  which  a  part  of  open 
elements  are  saturated).  Since  both  the  p  —  Vpop-open  a^ad  the  p  —  Vsat  relations  as  well  as 
p—V{p  >  PoVto)  equation  are  represented  by  a  common  function,  I3  (p)  (=  erf{Cp), 
C  =  ■v/i'A/4,  p  —  p/po  —  1))  the  pressure  at  the  inflection  point,  po,  is  not  only  the  pressure 
at  which  the  local  compliance,  dV/dp,  is  maximum,  but  also  the  pressure  location  for  a 
maximum  rate  of  increase  in  Vpop-open  as  well  as  in  Vsat  with  their  gradients  given  by  the 
following  equations; 

Po5to<P, 
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dVpop—open 

dp 


AVA 


1  +  Vto  4po 


exp{—{Cp)^)  for  p>0. 


(10) 


The  identity  between  dVsat/dp  and  dV/dp  implies  that  the  shape  of  the  p-V  curves  (in  the 
range  of  the  high  pressure  solution)  closely  represents  the  change  in  the  saturated  volume 
rather  than  the  pop-open  volume.  Since  the  gradient  of  Vpop-open  is  smaller  than  that  of 
Vsat  by  a  factor  of  1/(1  -f  yro)  (<  1),  the  magnitude  of  Vgat  eventually  becomes  greater 
than  that  of  Vpop-open  as  pressure  increases.  Also,  1/(1  +  yro))  a  fraction  of  the  pop-open 
volmne  (=  pop-open  volume/total  volume  of  a  single  TRS  element),  may  be  interpreted 
as  a  gradient  ratio  of  Vpop-open  to  Vgat-  A  smaller  value  of  yro  (i.e.  smaller  piston  stroke 
limit)  means  that  the  element,  once  it  pops  open,  reaches  the  saturated  state  earlier; 
hence,  for  Data  T  of  Fig.l2(f)  (yro  =  0.289,  PqVtq  =  8.68  [cmH20\)  Vsat  becomes  greater 
than  Vpop-open  at  a  pressure  close  to  po  yro,  while,  for  Data  M  of  Fig.l2(c)  (yro  =  0.626, 
PoVto  =  18.98  [cmH20])  Vsat  does  not  overtake  Vpop-open  within  the  measured  pressure 
range.  The  gradients,  dVpop-open/ dp  between  p  =  0  and  p  =  Pf  (=final  pressure)  (solid) 
and  dVsat/dp  between  p  =  po  yro  and  p  =  Pf  (broken)  are  plotted  in  Fig.13  for  the  six  data 
sets.  The  gradients  are  symmetric  with  respect  to  po-  The  data  sets  with  high  gradients 
(Data  set  M,  N,  T)  show  high  sensitivity  of  the  gradients  to  pressure  change  near  po-  Other 
data  sets  with  low  gradients,  particularly  Data  set  E  and  R,  indicate  that  the  gradients 
(i.e.  the  local  compliance)  do  not  change  too  much  over  a  substantial  range  in  pressure 
around  po- 
Summary 

The  mechanistic  model  of  TRS  developed  in  Part  I  is  applied  to  examine  p-V  curves 
(in  a  form  of  the  error  function  p-V  equation)  of  patients  with  ARDS  with  the  following 
results: 

1.  Parameters  of  the  deflation  process,  Vy  in  Eq.(4)  and  in  Eq.(5),  predicted  by 
the  mechanistic  model  of  the  corresponding  inflation  process  agree  well  with  those  of 
the  error  function  p-V  equation  for  the  deflation  process  (Figs. 2,  4),  indicating  that  the 
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mechanistic  model  has  a  certain  validity  to  be  used  for  improving  our  quantitative  un¬ 
derstanding  of  various  intra-respiratory  conditions,  that  the  shape  and  characteristics  of 
deflation  cirrves  depend  not  so  much  on  the  inflation  history  but  on  the  end-of-inflation 
(the  onset-of-deflation)  pressure,  and  that  relations  between  parameters  in  p-V  equation 
and  TRS  conditions  are  needed  to  strengthen  our  applications  of  p-V  curves  in  clinical 
settings  (Fig.  5). 

2.  The  non-dimensional  p-V  curve,  p  —  V  curve,  is  effective  in  distinguishing  differences  in 
magnitudes  of  model  parameters  among  different  p-V  curves  (Figs.7,  8). 

3.  In  the  mechanistic  model,  the  distribution  function,  Eq.(3),  and  its  change  with  pressure 
are  the  basis  for  evaluating  alveolar  recruitment  and  the  elastic  tissue  distension.  The  shape 
of  the  distribution  function  (the  peak  value  and  the  standard  deviation)  is  determined  by 
the  magnitude  of  the  non-dimensional  parameter.  A;  while,  the  magnitude  of  the  pressure 
at  the  ma.Yimnm  compliance,  po,  and  its  location  relative  to  the  range  of  the  p-V  curve 
are  the  important  factors  affecting  changes  of  the  distribution  with  pressure  (Eqs.(7)(9), 
Figs.9,10,11). 

4.  In  addition  to  A  and  po ,  other  parameters  of  the  model,  A V,  Vu  and  yro  influence  the 
magnitude  and  changes  of  both  Vpop-open  (volume  increase  due  to  alveolar  recruitment) 
and  Vsat  (total  volume  of  saturated  (fully-distended)  elements).  The  shape  of  the  clinically- 
measured  p-V  curve  represent  the  change  in  Vsat-  The  inflection  pressure,  po,  is  not  only  the 
pressure  at  which  the  local  compliance,  dV/dp,  is  maximum,  but  also  the  pressiure  location 
for  a  maximum  rate  of  increase  in  Vpop-open  as  well  as  in  Vsat  (Eqs.(8),(10),  Figs.11,12). 

The  mechanistic  model  of  a  TRS  element  presented  in  this  report  consists  of  a  sim¬ 
ple  piston-spring-cylinder  system  with  the  critical  pop-open  pressme  of  the  element  as 
distribution  parameter.  The  pop-open  volume  (=  volume  that  pops  open  at  the  critical 
pressure)  as  well  as  the  spring  constant  are  assumed  constant  and  common  to  all  elements. 
More  comprehensive  and  detailed  analyses  of  clinical  data  as  well  as  advice  from  clinical 
experts  are  needed  to  advance  the  model  further  and  also  to  make  it  a  practical  tool  for 
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understanding  various  respiratory  conditions.  However,  it  is  believed  that  the  analyses 
presented  here  show  the  developments  and  use  of  a  mechanistic  model  as  a  possible  new 
approach  to  investigate  respiratory  systems. 
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Figure/Table  Captions 
Table  1.  Summary  of  Deflation  Data  Analysis. 

Table  2.  Parameters  of  Inflation  Data  Sets  Examined 

Fig.  1.  Data  points  for  inflation  (unfilled)  and  deflation  (filled),  and  the  corresponding 
p-V  equation  of  mechanistic  model  for  inflation  (I),  and  the  error  function 
p- V  equation  for  deflation  (D) , 

PiD  —  pressure  at  the  intersect  of  the  inflation  curve  and  a  line  parallel  to  the  x-axis 
passing  through  the  initial  deflation  data  point. 

(See  Table  1  of  Part  I  (II)  for  numerical  values  of  parameters  for  inflation  (deflation).) 

Fig.  2.  Vy  (predicted  from  the  mechanistic  model)  vs  Vy  (of  error-function  p-V  equation 
for  deflation). 

Fig.  3.  VJy  (predicted  from  the  mechanistic  model  with  yro  =  0)  vs  (of  error-function 
p-V  equation  for  deflation). 

Fig.  4.  AV^  (predicted  from  the  mechanistic  model)  vs  AV'*  (of  error-function  p-V 
equation  for  deflation).  Letters  in  the  figure  =  Data  No. 

Fig.  5.  Deflation  cmrves. 

triangle  =  data  points,  dotted  =  error  function  p-V  equation  with  the  method  of 
least  squares  applied  to  determine  (A*^,  Pq,  AV*^,  Vjf), 

solid  =  error  function  p-V  equation  with  (AV'^,  V^)  determined  from  Eqs.(4,5) 
and  {A^,  Pq)  determined  by  the  method  of  least  squares. 

Fig.  6.  Ranges  of  parameters  of  inflation  data  sets. 

(a)  po  vs  A,  (b)  AV  vs  A,  (c)  1/(1  +  yro)  vs  A.  Letters  in  the  figure  =  Data  No. 

Fig.  7.  Model-based  p-V  equation. 

(a)  Data  Set  B,  E  and  N.  (b)  Data  Set  E  and  T. 
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(c)  Data  Set  N  and  R.  (d)Data  Set  E  and  M. 

Vertical  lines:  dotted  =  po>  broken  =  poyro- 

Fig.  8.  Non-dimensional  {p  —  V)  equation. 

(a)  Data  Set  B,  E  and  N.  (b)  Data  Set  E  and  T. 

(c)  Data  Set  N  and  R.  (d)Data  Set  E  and  M. 

Fig.  9.  Distribution  of  elements. 

(a)  through  (f)  for  Data  Sets  B(  =(a)),  E,  M,  N,  R  and  T(  =(f)). 

Fig.lO.  Niunber  fraction  of  total  open  elements,  Nopen  total  {p  =  P)/N,  vs  pressure  (solid) 
and  number  fraction  of  satmated  elements,  Ngat.  {p  =  p)/N,  vs  pressure  (broken). 

(a)  through  (f)  for  Data  Sets  B(  =(a)),  E,  M,  N,  R  and  T(  =(f)). 

Fig.ll.  Number  fraction  of  open  elements  at  p  =  0,  Nopen{p  —  0)/iV  vs  A 

Fig.l2.  p  —  V  curve  (solid,  Eq.(2)),  p  -  Vpop-open  curve  (dotted,  Eq.(8a)) 
and  p  —  Vaat  curve  (broken,  Eq.(8b))  for  six  data  sets. 

(a)  through  (f)  for  Data  Sets  B(  =(a)),  E,  M,  N,  R  and  T(  =(f)). 

Letters  in  the  figure  =  Data  No. 

Fig.l3.  dVpop-open/dp  [L/cmH20]  vs  p  [cmH20]  (solid) 

and  dVsatfdp  [L/cmH20]  vs  p  [cmH20]  (broken)  for  six  data  sets. 

(a)  through  (f)  for  Data  Sets  B(  =(a)),  E,  M,  N,  R  and  T(  =(f)). 


85 


Table  1.  Summary  of  Deflation  Data  Analysis 


Data  No 

Ad 

Po 

[cmH20] 

[L] 

Vi 

[L] 

vS 

[L] 

A  1. 

1.4712 

11.410 

1.9214 

-0.1528 

1.7685 

2. 

1.7567 

1.7064 

B  1. 

1.5528 

12.975 

1.4115 

-0.1227 

1.2887 

2. 

1.3405 

1.2808 

C  1. 

2.3293 

14.808 

1.1753 

0.0800 

1.2599 

2. 

1.2525 

1.2259 

D  1. 

1.0772 

6.846 

2.7012 

-0.6096 

2.0915 

2. 

2.5850 

2.0572 

E  1. 

2.2044 

15.602 

0.7009 

0.0703 

0.7712 

2. 

0.7784 

0.7329 

F  1. 

1.7995 

14.418 

1.0476 

-0.0527 

0.9948 

2. 

1.0041 

0.9252 

G  1. 

0.9839 

8.140 

1.2811 

0.0454 

1.3265 

2. 

1.6496 

1.3455 

H  1. 

0.7147 

6.472 

2.0772 

-1.0349 

1.0423 

2. 

1.1910 

1.0053 

I  1. 

2.1010 

14.157 

1.4798 

0.5637 

2.0435 

2. 

1.9943 

1.9784 

J  1. 

1.6600 

13.404 

1.8960 

-0.1250 

1.7103 

2. 

1.7759 

1.6858 

K  1. 

1.0577 

9.778 

1.5070 

-0.3657 

1.1412 

2. 

1.1620 

1.0860 

L  1. 

0.3107 

3.151 

1.3930 

-0.8076 

0.5853 

2. 

0.8574 

0.6206 

M  1. 

2.1702 

14.989 

1.4322 

0.1442 

1.5764 

2. 

1.4921 

1.5163 

N  1. 

0.7220 

6.45 

2.4020 

-1.1900 

1.2140 

2. 

1.1654 

1.1115 

0  1. 

1.0558 

10.969 

1.8684 

-0.3195 

1.5488 

2. 

1.6515 

1.4375 

P  1. 

1.5406 

12.231 

0.7461 

0.0684 

0.8145 

2. 

0.8679 

0.8154 

Q  1. 

1.1691 

10.670 

1.1801 

-0.1515 

1.0285 

2. 

1.0938 

0.9857 

R  1. 

0.8643 

8.1411 

1.8260 

-0.4893 

1.3366 

2. 

1.7784 

1.2923 

S  1. 

2.7458 

17.050 

2.0015 

-0.5673 

1.4341 

2. 

1.4461 

1.3962 

T  1. 

4.5056 

23.619 

1.1748 

0.1057 

1.2805 

2. 

1.2115 

1.2026 

U  L 

2.9735 

18.599 

2.3404 

0.0408 

2.3812 

2. 

2.5197 

2.3305 

1.  Obtained  by  method  of  least  squares  with  error  function  p-V  equation 
for  deflation  process. 

2.  Results  from  the  mechanistic  model  of  inflation  process. 
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Table  2.  Parameters  of  Inflation  Data  Sets  Examined 


Data  A 

Po 

AV 

1/(1  +  Vto) 

Vto 

<^D 

^open—sat.{P  — 

B 

2.7304 

21.999 

1.5559 

0.736 

0.359 

12.857 

0.01789 

E 

2.6497 

30.817 

1.6847 

0.676 

0.480 

18.559 

0.02645 

M 

2.9972 

30.327 

4.2463 

0.615 

0.626 

16.147 

0.04932 

N 

2.8046 

15.297 

1.6219 

0.736 

0.358 

8.704 

0.01685 

R 

1.6209 

16.352 

1.7396 

0.711 

0.406 

16.098 

0.07779 

T 

5.4708 

30.037 

1.7694 

0.766 

0.289 

8.761 

0.00012 

Po  and  (T£)  in  [cmH20],  AV  and  Vopen-satXP 

=  0)  in  [L]. 
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Fig.  1 .  Data  points  for  inflation  (unfilled)  and  deflation(filled),  and  the  corresponding  p-V 
equation  of  mechanistic  model  for  Inflation  (I)  and  the  error  function  p-V  equation  for  deflation  (D). 
PID  =  pressure  at  the  intersect  of  the  inflation  curve  and  a  line  parallel  to  the  x-axis  passing  through 
the  initial  deflation  data  point. 

(See  Table  1  of  Part  l(ll)  for  numerical  values  of  parameters  for  inflation  (deflation).) 


Vtq 
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(of  error  function  equation  for  deflati 


Fig  5.  triangle  =  data  points,  dotted  error  function  p-V  equation  with  the  method  of 
least  squares  applied  to  determine  {  ,Po,AV‘‘  ,V^) , 

Solid  =  error  function  p-V  equation  with  {AV‘‘,V^ )  determined  from  Eqs.  (4,5)  and  {A^,Pq) 
determined  by  method  of  least  squares. 
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H20] 


dNj/N  •  dpcj  [l/cmi?20] 


open  total  / ’iy  f 


p  [cmHtO]  p  [cmH^O] 


Fig.lO.  Number  fraction  of  total  open  elements,  (p  =  p)/N,  vs  pressure  (soUd) 

and  number  fraction  of  saturated  elements,  N,„.  (p  =  p)/iV,  vs  pressure  (broken). 

(a)  through  (f)  for  Data  Sets  B(  =(a)),  E,  M,  N,  R  and  T{  =(f)). 

9  7 
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Fig.ll.  Number  fraction  of  open  elements  at  p  =  0,  NopeniP  =  0)/N  vs  A  Fig .  1 1 
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Fig.l2.  p-V  curve  (solid,  Eq.(2)),  p  -  Vpop-open  curve  (dotted,  Eq.(8a)) 
and  p  —  Vsat  curve  (broken,  Eq^(8b))  for  six  data  sets. 

(a)  through  (f)  for  Data  Sets  B(  =(a)),  E,  M,  N,  R  and  T(  =(f)). 
Letters  in  the  figure  =  Data  No.  q 
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Fig.13.  dVj^j>-^en/dp  [L/cmH20\vs  p  [cmH20\  (solid) 

and  dVsat/dp  [L/cmH20]  vs  p  [cmH20]  (broken)  for  six  data  sets. 

(a)  through  (f)  for  Data  Sets  B(  =(a)),  E,  M,  N,  R  and  T(  =(f)). 
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Key  Research  Accomplishments 


1.  A  representation  of  p-V  ciirves  by  a  single  (non-linear)  model  equation  (either  the 
error-function  p-V  model  equation  or  the  sigmoidal  model  equation)  is  confirmed  to 
be  an  effective  method  for  clinical  data  analyses. 

2.  A  mechanistic  model  for  the  inflation  process  is  constructed;  which  makes  it  possible 
to  yield  information  on  the  internal  conditions  of  the  respiratory  system  firom  the  p-V 
curve. 

3.  The  ability  of  the  mechanistic  model  of  the  inflation  process  to  yield  information 
on  the  deflation  process  correctly  justifies  its  further  development  as  a  new  tool  for 
processing  p-V  curves  in  practical  environments. 
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Conclusions 


It  is  shown  that  both  the  sigmoidal  (tangent  hyperbolic)  p-V  model  equation  and  the 
error  function  model  equation  represent  quasi-static  p-V  curves  well.  Major  parameters  of 
both  the  sigmoidal  and  the  error  function  model  equations  are  the  non-dimensional  com¬ 
pliance,  A,  the  maximum  local  compliance,  po,  the  upper  (or  lower)  volume  asymptote, 
Vu  (Vl),  and  the  maximum  voliune  available  for  inflation,  AV.  Although  both  continuous 
model  equations  are  antisymmetric  with  respect  to  poi  the  non-dimensional  parameter, 
A,  as  well  as  two  volume  asymptotes  slightly  differ  between  the  two  equations  as  those 
function-specific  parameters  are  selected  to  follow  a  specified  i>-V  curve  as  closely  as  pos¬ 
sible.  The  p-V  model  equation  in  the  form  of  the  error-function  (because  of  its  relevance 
to  the  mechanistic  model)  is  applied  to  data  sources  of  patients  with  ARDS  as  well  as  of 
healthy  adults  to  show  that  the  the  magnitudes  of  the  model  equation  parameters  distin¬ 
guish  various  p-V  curves  clearly. 

A  mechanistic  model  of  TRS  elements,  each  consisting  of  a  piston-spring  system, 
is  developed  to  analyze  p-V  curves  for  the  inflation  process.  The  model  accommodates 
both  the  alveolar  recruitment  (in  terms  of  the  critical  pop-open  pressure)  and  the  elastic 
distension  of  wall  tissues  (in  terms  of  the  piston  displacement).  Model-based  relations  are 
established  between  the  parameters  in  the  p  —  V  curve  represented  by  the  error  function 
equation,  and  in  the  mechanistic  model.  The  p-V  equation  derived  from  the  mechanistic 
model  consists  of  two  equations;  one  for  the  low  pressure  region  where  all  open  elements 
are  active  (=  imsaturated)  as  the  piston  of  an  element  is  yet  to  reach  its  stroke  limit, 
and  the  other  for  the  high  pressure  region  where  some  open  elements  are  saturated.  The 
elemental  distribution  over  the  critical  pop-open  pressme  is  a  normal  distribution  with  its 
shape  (the  mean  and  the  standard  deviation)  affected  substantially  by  the  magnitudes  of 
two  parameters  of  the  model,  A  and  po- 

The  mechanistic  model  is  applied  to  p-V  curves  of  patients  with  ARDS.  Parameters 
of  the  deflation  process,  and  AV**,  predicted  by  the  mechanistic  model  of  the  corre- 


spending  inflation  process  agree  well  with  those  of  the  p-V  model  equation  for  the  deflation 
process  (Chapter  3,  Figs.2,  4),  thus  providing  justiflcations  for  its  validity  to  be  used  for 
improving  our  quantitative  understanding  of  various  intra-respiratory  conditions.  In  the 
mechanistic  model,  the  distribution  of  TRS  elements  and  its  change  with  pressure  are  the 
basis  for  evaluating  alveolar  recruitment  and  the  elastic  tissue  distension.  The  shape  of 
the  distribution  function  (the  peak  value  and  the  standard  deviation)  is  determined  by  the 
magnitude  of  the  non-dimensional  parameter,  A;  while,  the  magnitude  of  the  pressure  at 
the  mflYiTmim  compliance,  pQ,  and  its  location  relative  to  the  range  of  the  p-V  curve  are 
shown  to  be  important  factors  affecting  changes  of  the  distribution  with  pressure.  The 
other  parameters  of  the  model,  AV,  Vu  and  yro  influence  the  magnitude  and  changes  of 
both  Vpop-open  (volume  increase  due  to  alveolar  recruitment)  and  V^at  (total  volume  of 
saturated  (fully-distended)  elements). 

We  are  currently  in  the  process  of  developing  a  method  of  analyzing  the  deflation  curve 
based  on  our  mechanistic  model.  After  a  mechanistic  model  is  derived  for  the  deflation 
process,  it  is  possible  to  analyze  both  the  inflation  and  the  deflation  process  together;  which 
would  yield  a  greater  amovmt  of  quantitative  information  on  the  respiratory  conditions  than 
the  analyses  of  the  inflation  process  alone. 

More  comprehensive  and  detailed  analyses  of  clinical  data  as  well  as  advice  from 
clinical  experts  are  needed  to  refine  the  model  further  and  also  to  relate  predictions  by  the 
mechanistic  model  to  clinical  diagnoses.  However,  it  is  believed  that  the  analyses  presented 
here  show  the  developments  and  use  of  a  mechanistic  model  as  a  new,  effective  approach 
to  investigate  respiratory  systems  of  patients  with  acute  lung  injury. 
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